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ABSTRACT 


This study is part of our continuing program for the development of 
remote sensing for geologic applications. The main objective was to develop 
techniques and evaluate our ability to discriminate among surface geologic 
materials on the basis of their thermal properties as determined from HCMM 
data. Three test sites in the western U. S. were selected for study. The 
results of the invv'sstlgation are: 

1) Attempts to determine quantitatively accurate thennal inertia values 
from HCMM digital data met with only partial success due to the effects of 
sensor mlscalibrations, radiative transfer in 'he atmosphere, and varying 
ffleteorology and elevation across a scene. 

2) In most instances, apparent thermal inertia is an excellent qualita- 
tive representation of true thermal Inertia. For many applications this will 
reduce datn processing requirements with only minimal loss of Information. 

3) Computer processing of digital day and night HCMM data allowed 
construction of geologically useful images. At some test sites, more inform- 
ation was provided by Heat Capacity Mapping Mission (HCMM) data than Landsat 
data. Soil moisture effects and differences in spectrally dark materials were 
more effectively displayed using the thermal data. 

A) Future work will extend these results, and examine the multivariate 
use of HCMM data with Landsat and Seasat data. 
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I, INTRODUCTION 


The use of middle IR thermal data represents one aspect of our develop- 
ment of remote sensing techniques for geologic mapping, mineral exploration 
and mineral assessment. Considerable success has been achieved in discrim™ 
inating lithologic units, and in the separation of altered materials using 
various enhancement techniques with the reflected solar visible and near IR 
data from Landsat and aircraft multispectral scanners (Goetz et al. , 1975; 
Abrams et al. , 1977; Goetz and Rowan, 1981). Tliese results have been verified 
with field data obtained with the JPL Portable Field Reflectance Spectrometer 
(PFRS) and lab spectra of samples from the field. Multispectral middle infra- 
red data, acquired by aircraft scanners, have also been demonstrated to separ- 
ate some geological units, in particular igneous silicate rocks (Vincent, 

1975; Kahle and Rowan, 1980). However, many ambiguities in the remotely 
sensed data remain. The addition of thermal Inertia data to the multi- 
spectral analysis will make possible the resolution of some of these problems. 
Some materials with similar reflective and emissive spectral characteristics 
will be able to be differentiated by their thermal properties. Also, since 
thermal properties are representative of the top several centimeters of the 
surface material, it will be possible to sense the characteristics of the 
material through a thin weathered layer or stain. 

As part of this continuing program at the Jet Propulsion Laboratory, we 
developed a numarlcal model of the thermal behavior of surface materials 
(Kahle, 1977). Using aircraft data, we created the first thermal Inertia 
image of a part of the Earth's surface at Pisgah Crater, California (Gillespie 
and Kahle, 1977) , and demonstrated the utility of this image for geologic 
mapping (Kahle et al. , 1976). Our purpose with the present study has been to 
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apply this technique to Heat Capacity Mapping Mission (HCMM) and aircraft data 
to further develop and demonstrate the use of middle IR thermal data for 
geologic applications. 
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Figure 1. Full-frame HCMM visible image. Test sites are indicated by 
D (Death Valley), W (Walker Lane), and P (Pisgah Crater). 
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criminate among rock units whose compositional differences are obscured in the 
visible and near-infrared regions by thin surface coatings, and 5) distinguish 
between outcrops and alluvial areas and among alluvial deposits. 

The second area of interest is Death Valley, California (D), also shown 
in Figure 1. Within Death Valley and the mountain ranges that border it, a 
great variety of rocks representing virtually every geologic period are 
exposed. Most of the topographic features of Death Valley can be attributed 
to movement along n major fault zone, the Death Valley fault zone, which 
approximately coincides with the valley. The valley is the drainage sump both 
for surface and subsurface flow for at least 23,000 square km of desert. The 
fill in the valley consists mainly of coarse-grained, alluvial-fan deposits on 
the slopett, and fine-grained alluvium and evaporlte deposits at the lower 
elevations of the valley. 

Previous study of NOAA satellite VHRR thermal infrared scanner data from 
the Death Valley region had shown very interesting temperature patterns at 
night. The mountains and valley floor are cool with the Intermediate eleva- 
tions being warm. This temperature pattern, which varies seasonally, appeared 
to be related to the evaporation from the valley floor. 

The objectives of the study of this area were to; 1) develop a better 
understanding of the physics of the spatial and diurnal temperature variations 
of this region, 2) determine the utility of thermal data for ground water 
mapping in a desert region, and 3) examine lithologic discrimination in a 
complex area. 

The third study area shown in Figure 1 is the Pisgah Crater - Lavlc Lake 
region of California (P). Here, the objectives were to: 1) compare a thermal 

inertia map derived from HCMM data with an already existing thermal inertia 
map derived from aircraft data (Kahle, et al. , 1976; Gillespie and Kahle, 
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1977), 2) determine the repeatability of the process under different meteor- 
ologleal conditions, and 3) determine the effects of reduced resolution In 
separating major lithologic units* 

A fourth site, San Riifael Swell, Utah for which only one UCMM digital 
data set was obtained has been Included primarily for future comparison with 
an existing Landsat-Seasat study (Blom et al », 1981). 

Tl\e approach taken at the three primary sites was essentially the S'lme. 
Prior to the HCMM launch, day and pre-dawn aircraft scanner data were acquired 
at each site. Tliese data usually consisted of several bands of multlspectral 
visible and near-Tk data and a single broadband thermal channel. Before and 
during the flights meteorological data were acquired at the site by field 
parties. Hie meteorological data were used in the JPL thermal model to derive 
tables of thermal Inertia as a function of day-niglvt temperature difference 
(and sometimes topography). Hiese tables were used with the aircraft scanner 
data to derive thermal Inertia images. 

tXirlng thf useful lifetime of HGMM, a similar procedure was followed. 
Picld measurements of meteorological data were made before and during known 
overpass times of the satellite. Tliese data were augmented by data from the 
closest National Weather Service stations. Using the thermal model and the 
meteorological measurements, the satellite data in digital format were reduced 
to thermal Inertia Images. 

Tliese thernuil inertia images, and also the temperature and albedo Images, 
were then evaluated for geological Information. Additional field verification 
needs to be undertaken. Effects of varying meteorological conditions, cali- 
bration problems, topographic corrections (both aspect and elevation), and 
atmospheric attenuation and radiance, have been evaluated. Tlve utility of the 
apparent thermal Inertia as compared with the more quantitatively correct 
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derived tltermal Inettia has been studied. Ihe inethods, results and Interpret- 
ation are discussed In the subsequent sections. 
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nX. DATA ACQUISITION 


Data were coilected from satellite and aircraft flights over the three 
test areas: Walker Lane, Nevada; Death Valley, California; and Plsgah Crater, 

California. Tlie three areas are shown on Images obtained from HCMM flights in 
Figures I through 4. The first Inuige (Figure 1) contains all three areas on a 
single, full-frame picture of the western United States. Hie subsequent 
figures show individual areas. In the Walker Lane area we have studied two 
test sites, the Goldfield mining district (G) and the Cuprite mining district 
(C) near Stonewall Mountain; both sites are indicated in Figure 2. Death 
Valley is showji in Figure 3 a’^d Plsgah in Figure 4. 

A . Aircraft Data 

Table 1 lists the digital aircraft data collected for this stiidy. 

Tlie aircraft data consist of nineteen different flights. Hie date of each 
flight, the area of coverage, the aircraft type and sensor used are given in 
the table. Double coverage of the Goldfield area was obtained on August 8, 
1977 with data from both the NP-3A and the U-2. Numbers in parentheses indi- 
cate more than one flight over an area. 

The aircraft operations consist of two flights over a test area 
within a 24 hour period. The two consecutive flights are necessary in order 
to determine the day-night temperature difference from the IR radiance data. 

A predawn flight is made Just before sunrise to measure the minimum diurnal 
temperature of the ground and a mid-afternoon flight, one or two hours after 
local solar noon, is made to measure the temperature maximum and the surface 
albedo. On some flights the NP-3A and NC-130B made two runs for each flight: 
one at a low level (~ 2.1 - 4.2 km ASL) and one at a high level (~ 7.6 km 
ASL) . The lengths of the flight lines are 3.7 km at Plsgah, 15-22 km at Death 


8 


ORIGINAL' PAGP 

BLACK AND WHITE PHOrOGRAPH 



Fl^ui^ 2. 


HCMM visible imr\yo of Walker Lane, Nevada. lest sites indicated 
by G (Goldfield mininH district) and C (Cuprite mining district). 
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Figure 3. HCMM visible image of Death Valley, California 
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Figure A. HCMM visible Image of I’lsgah Crnlcr, California. 
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TABtE I 


AIRCRAFT FLIGHTS 


Date 

Site 

Aircraft 

Sensor 

March 29-30, 1975 

Pisgah (2) 

NP-3A 

M^S 

September 3, 1975 

Pisgah 

U-2 

HCM 

June 14, 1976 

Walker Lane 
Death Valley 
Pisgah 

u-2 

HCH 

October 26-27, 1976 

Goldfield 

Cuprite 

NP-3A 

M^S 

March 11, 1977 

Death Valley 

NP-3A 

M^S 

July 13-14, 1977 

Pisgah 

NP-3A 

M^S 

August 8-9, 1977 

Walker lane 

Goldfield 

Cuprite 

U-2 

NP-3A 

HCM 

M'^S 

February 2-3, 1978 

Death Valley 

U-2 

HCM 

March 14-15, 1978 

Death Valley 

U-2 

HCM 

August 18, 1978 

Goldfield 

Cuprite 

NC-130B 

M^S 

April 4, 1979 

Death Valley (2) 

U-2 

DMS 
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Valley, and 5.5 km (fot the Goldfield area) and 8.5 km (foe the Captite area) 
at Walker lane. All f limits were at 19.8 km ASl. 'Ihe Bendlx M^S (Modular 

Multiapectral Scanner) was uaed on the Nl'-3A and the NC'"130B. This Is an 
electro-optical scanner which measures electromagnetic radiation In ten 
channels between 0.33 and l.l Mm, and In a thermal IR channel from 7.95 to 
13.5 Mm. The Instrument has a scan angle of 100^, a scan rate of 20 revolu- 
tions per second, and an Instantaneous field of view (IFOV) of 2.5 mllll“ 
radians. At an altitude of 2.1 km, the sensor covers a ground swath of 5.0 km 
In width and at nadir has an IFOV of about 5 m squaiA. The U-2 was flown with 
either of two sensors, the Heat Capacity Mapper (HCM) or the Daedalus Multl- 
spectrai Scanner (DMS). The HCM Is a scanner that was built at the Goddard 
Space Flight Center (G3FC). It detects the visible and near IR (0.5 — 0*7 Un) 

and thermal IR (10.5 - 12.5 Mm) portions of the electromagnetic spectrum. The 

HCM has a 90° scan angle, a 2.8 mllllradlan IFOV, and a thermal resolution of 

0.2°C. At an altitude of 19.8 km, the sensor scans a swath of 39.6 km, and at 
nadir has an IFOV of about 55 m square. The DMS Is an eleven channel scanner 
with a scan rate of 10 revolutions per second which covers the spectrum from 
0.38 to 1.1 Mm and from 10.4 to 12.5 Mm. The system uses either of two scan- 
heads of 1.25 or 2.5 mlllltadians angular resolution. From an altitude of 
19.8 km, the sensor covers a swath of either 15 or 36 km and has a resolution 
of either 25 or 50 m depending upon the scanhead. 

B. Satellite Data 

The HCMM satellite views the Earth's surface with the HCMR scanning/ 
imaging radiometer which operates with two channels in the visible and near IR 
(0.55 - 1.1 Mm) and the thermal IR (10.5 - 12.5 Mm) regions. The sensor has a 
scan rate of 14 revolutions per second, an IFOV of 0.83 mllllradlans, and a 
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thermal resolution of 0.4*^K at 280°K. At the satellite altitude of 620 km, 
the ground resolution is about 500 m. 

Prior to what appeared might be favorable weather conditions during 
satellite passes over our test sites, a team of Investigators was dispatched 
from JPt to obtain ground truth measurements; however, these measurements were 
only taken at one of the sites on any given occasion. These measurements are 
described in detail in the next section. Hourly weather observations from 
nearby National Weather Service stations were used as ground truth for the 
other two sites. Ground truth measurement"' at the site continued only if the 
weather remained clear; otherwise, the effort was abandoned. All the HCMM 
data of the test sites were first obtained from GSFC in image form and then 
examined to select only those passes when weather conditions permitted clear 
views of the ground. Digital HCMM images were obtained for periods when field 
measurements had been made and/or when the GSFC Images Indicated that good 
quality satellite data were available. Ilie digital HCMM data used in this 
study are listed in Table 2. 

HCMM digital data were obtained for 9 sets of overpasses from May, 
1978 to April, 1979 for the Walker Lane, Death Valley, and Plsgah Crater 
areas. With one exception (July 21-22, 1978), each of these sets includes a 
night overpass (N), 24 hours later another night overpass, 12 hours later a 
day overpass (D), and 24 hours later another day overpass. Each NNDD set 
covers most of all three test sites. Included in this data set is one com- 
pletely cloud- free nine-day period between July 15 and July 23, 1978, when 
field measurements were taken at Cuprite for the entire period and two HCMM 
NNDD sets of overpasses (7 total) were obtained. Simultaneous HCMM and air- 
craft data are available for April, 1979 flights over Death Valley. 
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TABT^K 2 


satkijU'i’K flights 

l)fl k c* H ^ 

May 13- IS, 1978 
May 29-31, 1978 
July S-7, 1978 
July 16-18, 1978 
July 21-22, 1978 
Auf»,ust 17-19, 1978 
AuHUHfc 27-28, 1978 
Scpkomhor 18-20, 1978 

Api'U 3-S, 1979 
C. .GtN ^um l 

A ground truth program wan oatablialu'd to oollcct supportlaig data 
for aircraft and HCMM overflights. The data consist of; 1) surface temper- 
ature measurements, 2) suhsurface temperatvire nHiasurements, 3) soil moisture 
measurements, and A) meteorological observations. The surface and subsurface 
measurementH were collected at various locations at each site, by a team of 
Investigators, and recorded by hand. The meteorological data were obtained in 
digital and analog formtit by autoiTMted eqnlpmonk placed at a central location 
at each site. For convenience, we. first discuss the surface and subsurface 
measurements - data sets (i) through (3), and then the meteorological 
measurements - data set (A). 


Remarks ^ ^ 

Vlsgah missed on 5/ I A night 

Death Valley and Flsgah missed on 3/30 night 

Flsgah missed on 7/6 night 

Death Valley and Flsgah missed on 7/17 night 

No second day coverage (duly 23) 

Plsgali missed on 8/18 night 

Han Rafael coverage only 

No Death Valley coverage 
Flsgah missed on 9/19 night 

All sites covered 
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Surfivco and BubHutfacc mnagnrninenfeH 


I1\e Hurfacn tempnvatnrn flieaguremonts wove* obfcalund wtbb a 
Barnofi VR'l-5 Precision foidiation 'nicrmomntnr* 'Ovls tnfitrument Is a portable 
infrared radiomett't whid\ operates In Cbe thermal region from 8 to 14 baa 
a field of view of about 2^^ and Is accurate to within 0.5^'C over a temper'- 
aturc( range from -20 to 175^0. 'ri\e radiation emitted from the target la read 
di rectly from a Rcale on the PR'f-S as an equivalent blackbody temperature. 

Subsurface tomperature measturements were made with Yellow 
Springs luHtrvuneni (Y»Sl) thermi.stora buried Ln the au,li at: a number of depths, 
usually 2 to 50 cm below the surface. These temperature sensors consist of 
vinyl-;|acketed probes whose tips contain thermistors which can be used over a 
range from -30 to HOO^C* Tlve probes are connected to a phone jack by leads, 
and when the jack is inserted in t)>e roceptacie of a YSl meter, the temper- 
ature is read from a scale. Figure 5 shows a diurnal pattern for soil temper- 
ature that was obtained in lUirch, 1978 from YSl probe measurements at IXiath 
Valley on a silty surface. Note that tbs surface heating maximum leads the 
maximum at the lower depths (- 25 cm) by 6 hours, and that the damping depth 
of the diurnal heating is about 50 cm. Both features are fairly typical of 
most soils. 

Tlie soil moisture measurements were obtained by the gravimetric 
method as follows. Samples are removed from the soil by an extractor, usually 
a cylindrical coring device. Ti\e moist soil is placed in a container, sealed, 
and later dried in an oven at about 60®C. Tl\e sample is weighed before and 
after drying to determine the moisture content. Uve soil moisture content can 
then be expressed in either of two ways* 1) as c percentage of the dry 
weight, (wet weight - dry weight) /dry weight, or 2) as a percentage of the 
volume, t molstvire by weight multiplied by the bulk density of: soil (dry 
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welght/field volume of soil sample). Figure 6 shows an example of volumetric 
soil moisture measurements for different depths on the floor of Death Valley 
opposite Desolation Canyon. Ihe ordinate is a scale indicating per cent of 
moisture content by volume. The abscissa denotes site locations starting at 
(10) which is furthest from the center of the Valley floor and ending at (6) 
which is nearest. Sites (10) and (9) are on the gravel fan beneath the 
canyon, (5) is on the silt past the fan, and (6) is on the salt pan which 
covers most of the valley floor. Sites (8) and (7) are located on transition 
zones between the silt and the fan, and the silt and the salt, respectively. 
The diagram shows the dryness of the fans, and the increase of moisture with 
distance along the silty surface until the salt pan is reached. 

2. Meteorological Measurements 

Meteorological measurements are needed to specify upper (or 
surface energy) boundary conditions for the soil thermal model. As an ex- 
ample, Figure 7 shows the surface energy balance for a day in July, 1978 at 
Stonewall playa which is located near Goldfield, Nevada. These are typical 
diurnal patterns in solar radiation (S), soil heat conduction (G), net thermal 
(longwave) radiation (R), and sensible heat flux (H) for a dry, desert 
surface. Hie methods used to calculate these values from periodic near- 
surface measurements are discussed in Kahle (1977). 

Radiation measurements, S and R, consist of net shortwave, and 
longwave flux densities. Net radiation flux densities are obtained by use of 
a Thornthwaite Model 602 radiometer and strip chart recorder. The radiometer 
is a flat, circular, black plate (thermal transducer) that is shielded by two 
clear plastic hemispheres, and has a detection range of about -280 to 

r\ 

flllO w/ra'^ where positive refers to a net radiant flux directed toward the 
ground and negative vice-versa. The shortwave component (0.28 - 2.80 jm) is 


18 


DEATH VALLEY - DESOLATION CANYON FAN TRAVERSE 


ORIGINAL PAGE IS 
OF POOR QUALITY 



^ ^ ^ a s ^ ° 
(%)iNaiNOD aynislow Diyi3wmoA 


.19 


Figure 6. Volumetric soil moisture variation with depth at Death Valley on March 15, 19/8* 




Figure 7. Surface heat flux density at Stonewall Playa, near Goldfield, Nevada 
for a day during July, 1978. 






measured with an Kpplcy Precision Spectral Pyranometer (PSF) and the. longwave 
component (S-SO pm) with an Bpplcy Precision Infrared Radiometer (PIR), All 
radiometers are mounted at a height of 1/2 to 1 meter above the ground, and 
the PSP or PIR can be used to measure either upward or downward radiative flux 
densities. 

Soil heat flux densities (G) are measured with a Soil Heat Plux 
Recording System, Model 310. The sensor is a Thornthwaite thermopile enclosed 
in a thin, flat disk - 25 mm in diameter and 2.6 mm thick. Tlic disk is placed 
just below the soil surface and parallel to it and the temperature difference 
across its vipper and lower faces is proportional to the soil heat flux. The 
data are entered on a strip chart recorder. 

Sensible (H) and latent (L) heat flux densities are computed 
from weather observations of wind speed, air temperature, and air humidity. 

In situ weather data are obtained at each site from two systems. Hie first 
system is a Meteorology Research, Inc. (MRI) mechanical weather station. This 
is a self-contained instrument powered by four 1.5 V flashlight batteries, and 
is equipped to measure wind speed and direction, air temperature, and relative 
humidity at a single level above the ground (~ 1 or 2 ra) . Data are recorded 
on a strip chart. The second system is a ralcroraeteorological station - 
designed and built at JPL - that collects data at 7 levels above the ground 
(1/8, 1/A, 1/2, 1, 2, A, 8 ra). A psychroraeter-aTiamometer pair are mounted at 
each level at the end of a 1 m rod that is attached to a mast a little over 
8 m tall (Figure 8). The anemometers are three-cup type about 9 cm high, 5 cm 
In diameter, and with rotor arms 7.5 cm long. The rotation rate, in RPM, of 
the cup assembly is related to the. wind speed (U) by the formula 
U * Q.026 x RPM f 0.118 rn/s which was obtained by wind tunnel tests. The 
relation between U and RPM - obtained by fitting a least squares line through 
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the data - is almost perfectly linear, r^ * 0.992 (where r is the correlation 
coefficient and “ 1.000 denotes perfect linearity). Tl>e psychrometers are 
a dry-bulb and wet-bulb arrangement of paired thermistors, Yellow Springs 
Instrument Model //44202, which are accurate to within 10.15®G over a range of 
-5 to 145'^C. Hie wef-bulb thermometer is covered by a cotton wick that is 
supplied with water from a small reservoir attached to the bottom of the 
psychrometer case. A small fan set in front of the case ventilates the therm- 
istor couple. The entire system is powered by a 6 volt battery. Data are 
printed In digital form on paper tapes, and can be averaged over several time 
intervals from 30 seconds to 16 minutes. 

D . Thermal Inertia Meter 

The paucity of laboratory measurements of the thermal inertia com- 
ponents (conductivity, density, and specific heat) for rocks and soils found 
in natural field conditions is apparent after review of the heat conduction 
literature. The vast majority of these measurements have been of engineering 
or insulating materials and not geologic materials in the field. For example, 
the Handbook of Physical Constants (Clark, 1966), presents approximately 75 
rock type conductivity values, only a few density values, and no specific heat 
values . 

Several methods of laboratory or direct in situ measurement of 
thermal Inertia have been developed. For example, Brennan (1971) has develop- 
ed a laboratory method of detorrralning thermal diffusivity (k*K/pc). It 
utilizes an Instrument having a transltor heated base plate which remains at a 
constant temperature through thermostatic control. A small rock core is 
placed on the plate and s thermistor monitors the temperature at the opposite 
end of the core. Thermal diffusivity is inversely proportional to the time 
taken by the unheated end of the core to rise an increment of temperature. 
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Bocnuso the 1 nstrumont ts calibrated, the diffusivity can then be calculnted 
when the temperature increment Is smjiU with respect to the teiriperature 
difference between the base plate and the initial core temperature. Brennan 
(1971) has published approximately forty measurements of thernril dlffusivity 
and thermftX inertia. 

A laboratory method has been developed by SchuXt?; (1968) for non^" 
destructively determining the. tbcrnuil inertia of materials at ambient temper'- 
ature.8. 1110 . method involves the radiant heating of an unknown sample and a 
standard. 'Ihe temperature history of the sample, and standard are monitored 
with an Infrared radiometer and comparison of their temperature histories 
yields the local thernkal inertia of the sample. 

In this study, tVte surface temperature rlvsc of the materials is used 
to calculate the thermal Inertia of the sample. Tills calculation is based 
upon the following assunptionst X) the compositions of the, standard and 
sample materials are homogeneous, 2) the heat flux density at the surface of 

each tnaterlal is constant and equal for each roateriai., and 3) the laaterlals 

arc thick enough to be approxlnuited as semi -inf Inite bodies. Tliis allows us 
to calculate the thernwl inertia of the sample from a simple, algebraic 

equation adapted by Schultz; from Carslaw and Jaeger (1959), 




‘^’'s 2 

* ’’s * • 

t 

where. ; 


^t 

target thermal Inertia 


- 

standard thermal inertia 


AT^ - 

standard change in temperature 


fT 

It 

target change in temperature 


^t " 

target emlssivity 



standard emisslvlty 
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A crlticfil point Is whetbet the net flux density at the surface Is 

constant# Even under laboratory conditions, the actual surface boundary 

condition Is non'-llnear because the surface is emitting radiant energy at a 

rate proportional to the fourth power of the surface temperature, 

Under controlled circumstances, l.e. in the laboratory, the surroundings 

radiate energy onto the test surfaces by an amount proportional to the fourth 

power of the room (ambient) temperature, so that the net energy loss 

4 4 

from the surface depends upon the difference, ' *^ 3 nib* '^amb 

usually changes slowly,, It can be cone'ldered constant and the non-linearity of 
the boundary condition depends almost entirely on the com|)onent, 

Provided the magnltuQc of the constant flux from the heat source is large and 
the duration of heating Is not too long, the net flux density at the surface 
will be approximately constant and the heating history of the materials will 
be represented by a parabolic curve (cf. Carslaw and Jaeger, 1959, pg. 75). 

wO 1 

Schultz has achieved these results by heating materials at 6 cal cm s” for 
5 seconds. My attempt to employ this technique In the field must deal with a 
more complicated set of environmental conditions, viz. the most serious dif- 
ficulty being gusty winds which transfer heat to, or from, the materials by • 
advection or convection. 

Nonetheless, the fundamental simplicity of Schultz's technique 
suggests that this method can be applied to construct a field Instrument to 
measure thennal inertia In situ . The JPL thermal Inertia meter (TIM) is 
designed to heat a target and two standards, with the resulting temperature 
change of each material monitored with an infrared radiometer. This permits 


' TIM was designed and tested by Dr. Stuart Marsh during his tenure at JPL 
as a National Research Council Research Associate. 
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the relationship expressed by Equation I to be used under field conditions. 
Again^ this relationship is valid only if the three assumptions listed above 
are satisfied, 'Hie first two conditions can be satisfied appraximately by 
proper selection of test materials » by calibration of the heat sources, and by 
monitoring the local environmentaJ conditions-prlncipally the wind gusts, 'fhe 
third condition is satisfied by insuring that the thickness of the standard 
material and target is greater than the damping depth of the heating wave. 

A photograph of the thermal Inertia meter and its configuration Is 
shown in Figure 9. The radiant energy sources (A) are 3200 W quart?, infrared 
heat lamps (OE #3200T3/1CL/HT-384V) . The lamps are njounted in parabolic, 
polished aluminum reflectors coated with lacquer to prevent disct'iloratlon. 

Each lamp is pivoted which permits it to be quickly moved into position to 
heat the surface and removed to allow monitoring of the surface temperature. 
The lamps are connected to Variac voltage regulators (B) which permit the 
mutual balancing of the three Infrared lamps. The surface temperature rise is 
monitored using a Barnes Engineering Company Precision Radiation Thermometer 
(C), model PRT-5. The Instrument is specified to be accurate to 0.5°C. The 
instrument's detector head (D) is mounted in a metal block in such a way as to 
permit rapid movement from reference standard to target. The detector head is 
38.7 cm above the target or standard yielding a circular field of view 1.35 cm 
in diameter at the center of the area undergoing heating. A portable digital 
voltmeter (Danameter)(E) measures the voltage from the PRT-5. The PRT-5 and 
voltmeter are mutually calibrated against a temperature-controlled blackbody. 

The base of the thermal inertia meter consists of three compart- 
ments. Their interiors are painted with flat-black, and are surrounded by 5 
cm of styrofoam insulation. Compartments 1 and 3 hold the thermal inertia 
standards. These standards are placed within plywood boxes, also painted with 
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Figure 9. Thermal Inertia Meter (TIM) with quartz infrared heat lamps (A), voltage regulators (B) , 
PRT-5 (C), and detector head (D) , digital voltmeter (E) , compartments for thermal inertia 
standards (1 and 3), compartment open to the surface to be measured (2). 



flat-black, and then positioned into the compartmenta. Compartment 2 is left 
open to the surface which is to be measured by the TIM. 

A dolomite block (30 x 18 x 5 cm) and 20/30 mesh Ottawa sand 
(quartz) were acquired for use as thermal inertia reference standards. These 
materials were selected in order to have both high and low thermal inertia 
references. Both materials were sent to an independent laboratory^ for deter- 
mination of their respective thermal conductivity, specific heat, and density. 
The specific heat of both standards was determined by the method of mixtures, 
the density by a Beckman air pycnometer, and the conductivity by the American 
Society of Testing Materials (ASTM) method. Results of these determinations 
are given in Table 3. 

Laboratory testing insured that each lamp was heating evenly across 
the surface of the standard or target. PRT-5 measurements indicated that the 
surfaces were being evenly heated within the accuracy of the PRT-5 (0.5°C). 

The V ^diative energy outputs of the three infrared lamps were balanced by 
adjusting the settings of the voltage regulators so that after a four minute 
heating period the temperatures of a given standard recorded during cooling 
differed by no more than 0.5°C for each lamp. This calibration procedure was 

TABLE 3 

THERMAL INERTIA STANDARDS 



DOLOMITE 

OTTAWA SAND 

Conductivity (R) cal nm”^s“^°C“^ 

0.01182 

0.00081 

Density (p) g cm~^ 

2.97 

1.74 

Specific Heat (c) cal 

0.20 

0.25 

Thermal inertia (P) cal cm~^ 



S-1/20C-1 

0.0838 

0.0188 


^(Geosciences LTD., Solana Beach, California 92075) 
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porf'ormGd with both the Ottawa sand and the dolomite standards to insure 
accuracy. 

Solution of Kquation 1 requires a hnowledRe of the infrared 
emissivities of the standards and target. Hie range of infrared emissivities 
found among natural terrain nviterials is generally between 0.89 and 0.99. If 
no correction is made for emissivity of target and standard in Kquation I, 
given a poasible range of 0.89 “ 0.99, an error of 35% could arise in the 
value, of tlietmal inertia. Wo. did not have the capability to measure the emis- 
sivity of the dolomite and Ottawa sand and on a practical basis it would be 
impossible to measure the emissivity of each target in the field. However, by 
employing published values (Buettner and Kern, 1965) for dolomite (0.96) and 
quartz sand (0-92), and by assuming a median value of 0.92 for the target, the 
potential error in calculated thermal Inertia is less than about 14%. 

Further testing of the TIM under laboratory conditions disclosed a 

number of problems. Tlie chief difficulty was determining the nature and 

magnitude of the heat flux density (Fq) at the surface of the standards and 

target. Use of Kquation 1 in the ana.lysls required that Fq be constant and 

the same for each material (cf. assiumption 2 in the previous discussion). 

Attempts to measure Fq with the Eppley precision IR radiometer under 

controlled, laboratory conditions proved difficult. A warm-up run of 

approximivtely seven minutes was necessary before Fq attained a constant value, 
-'2 —1 

about 0.022 cal cm ^ for each lamp. 

Tliis value of Fq was comiiared with numerical calculations of F 
derived from measurements of the heating and cooling rates of dolomite and 
Ottawa sand based on detailed laboratory studies. TiTe surface temperature 
histories of the two materials were calculated from a finite difference 
approximation to the one-dimensional heat conduction equation. Next, the 
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surface boundary condition was adj^isted until the calculated surface 

temperature at the end of the heating phase matched the observed temperature 

at this point, and then the surface temperature was calculated for a three 

minute period of cooling following the removal of the heat source. Calculated 

values of the cooling history were compared with the observed temperature 

decreases (Figures 10 and 11) in order to infer the nature of the surface 

boundary condition. In these calculations, the surface boundary condition was 

comprised of two parts, a constant flux Fq and a term accounting for radiative 

heat loss at the surface, o (T^ “ T^), where o is the Stefan-Boltzmann 

constant. Other trials were made with only the Fq term included in the 

calculations, but the results were no better than those shown here. Even so, 

the calculated values of Fq in Figures 10 and 11 were both larger (0.025 and 

0.030 cal cm“^s”^ respectively) than the measured value of 0.022 cal cm ^s”^, 

and furthermore, the calculated values were not equal. At first, it was 

thought that some heat loss due to convection could be occurring at the 

surface, but incorporation of an empirical ex|)resslon (Hampton, 1946) to 

account for this effect did not appreciably alter the results. The calculated 

results for the heating phase were even more inconsistent than those for the 

cooling phase. This is shown in Figure 12 in which the observed increase in 

surface temperature for Ottawa sand is compared with the calculated increase 

for three different boundary conditions. Values of Fq ranged from 0.042 to 

0.057 cal cm s , about a factor of 2 more than the measured value. A 

similar calculation for dolomite produced an even larger value of Fq, 0.070 
-2 -1 

cal cm s . 

Despite this difficulty, there is some merit in using the TIM, in the 
field to distinguish different materials from each other based on relative 
differences in the estimated thermal inertia. As long as wind speeds are low. 
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Figure 10, Model calculation for Ottawa sand. 1 'q is choseii to predict the 
correct temperature at the end of the heating period. 

31 



TEMPERATURE ( C) 






PAGE IS 
OF POOR QUALITY 



Figure 12. Calculated heating curves for Ottawa sand compared with observed 
heating after a warm-up period of ten minutes. 


33 



80 that convective heat .loss from the surface is mlnimli'-ed, the assumption of 
a constant flux at the surface under each lamp may be warranted. To test this 
hypothesis, we have used the TIM on field trips in the western United States. 

Measurements were obtained during two field trips to selected test 
sites in western Nevada and eastern California in August and December 1979. 

The field sites in and around Goldfield, Hawthorne, and Rlko, Nevada, and 
Flsgah Crater, California, were chosen because thermal modeling or thermal 
inertia mapping either had been performed or was planned over these sites. 
IXirlng field trips to these sites, approximately 60 measurements with the 
thermal inertia meter were recorded. As thermal inertia is very sensitive to 
soil moisture variations (Qulel, 1975), soil moisture samples of all alluvial 
material were taken. All sites were located in arid cliim»tes and the soils 
were quite dry (less than 10% moisture content per nnit volume near the sur- 
face). Of these 60 measurements, approximately 20% were unusable due to 
equipment malfunction or strong gusts of wind during the heating or cooling 
cycle. 

Table 4 shows the results of this investigation. Tlie value of the 
calculated thermal inertia and the error are given for several rock and soil 
types based upon both the Ottawa sand and dolomite measurements# The error 
limits are only meant to Indicate the effect of the 0.5°C inaccuracy in the 
PRT-5 measurement. A previously published value for a similar material is 
also given when available. Note that the dolomite-based values are consist- 
ently larger (about 21% on average) than the Ottawa sand-based values. Part 
of this discrepancy is due to the different heating histories of both 
materials. Both samples and targets were heated for less than seven ral.nutes 
so that the heat flux density at the surface was still changing. Another 
factor is that we used typical values of emlsslvlty for both standards. Any 
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TABLE 4 


thermal inertia derived from in situ thermal 


INERTIA METER MEASUREMENTS 


Material 

P Relative to 
Ottawa Sand 

P Relative to 
Dolomite 

Published 

Value 

Olivine Basalt 

0.032 

Jt 

0.003 

0.042 

k 

0.005 

0.053 

(Janza) 

Pahoehoe Basalt 

0.027 

k 

0.002 

0.039 

k 

0.004 



Basalt 

0.032 

k 

0.002 

0.042 

k 

0.003 



Barite 

0.038 

k 

0.003 

0.043 

k 

0.005 



Cliert 

0.042 

k 

0.005 

0.053 

k 

0.005 



Andesite 

0.033 

k 

0.003 

0.044 

k 

0.0C3 

0.055 

(Brennan) 

Rhyolite Ash Flow Tuff 

0.016 

k 

0.001 

0.022 

k 

0.003 



Sllicified Rhyodacite 

0.040 

k 

0.004 

0.048 

k 

0.005 



Aeolian Sand 

0.012 

k 

0.001 

0.015 

k 

0.002 



Clay -Si It Playa 

0.018 

k 

0.001 

0.024 

k 

0.002 



Sandy Alluvium 

0.010 

k 

0.001 

0.014 

k 

0.002 

0.014 

(Clark) 


deviation from the actual emtssivlties would cause a consistent difference 
between the P value based on the dolomite standard and that based on the 
Ottawa sand standard. For example, an error of 5% in either standard will 
introduce a consistent error of 10% between the two calculations. Nonethe- 
less, the table shows that it is possible to discriminate materials like sands 
and silts from barites and basalts. Also, the ratio of P values for olivine 
basalt and andesite, versus sandy alluvium is nearly the same (~ 3.5) for both 
published and calculated values. 

The TIM field tests indicate that the instrument is useful in a 
relative sense. As long as very accurate measurements of thermal Inertia are 


35 




not required, it is adequate i'/r discriminating among different materials. 

For example, our measurements of thermal inertia near Rlko were to test the 
premise that exploration for barite might be feasible because of barite's very 
high density and hence high thermal inertia. We used the TIM to measure the 
thermal Inertia of some very high grade barite deposits and the surrounding 
rock materials. Somewhat surprisingly the chert, which was commonly in the 
same area as the barite, in all cases had a higher thermal Inertia than the 
barite. This was true both for outcrops and rubble. Wliile barite has a 
higher density than chert, chert must have a higher thermal conductivity. 
Regardless of any error in the absolute value of thermal Inertia of these two 
materials, the use of the TIM demonstrated that the chert consistently had a 
higher thermal inertia than the barite. Both had higher thermal inertias than 
most of the surrounding rocks, so thermal inertia Images could probably be 
used to trace the extent of the chert-barite units. 

Tlie TIM is not useful in discriminating between material?? P 

values are nearly the same nor for absolute measurements of thermal inertia, 
yet the need for such a field instrument is recognized. For this use, the 
precision of the instrument must be enhanced. Further study and testing is 
required, particularly in the following areas: 1) rapid application (~ 1 min 

or less) of a large, uniform and constant heat flux to the standard and target 
surfaces, 2) a detailed study of the best kinds of materials for use as 
standards, and 3) the effect of convective heat loss in the calculations. 
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IV. DATA PROCfiSSINO 


Tlie reraotci sensing and gcovind truth data are used in the JPL thernutl, 
model as elements of a digital algorithm to compute thermal inertia values. 

The model is a one-dimenRlonal heat conduction equation which is solved by 
finite difference techniques, subject to a specific set of initial and 
boundary conditions. In its original form, the model describes changes in 
temperature with depth and time for dry soils - an application that is 
reasonable for arid climates. A second version of the model has been 
developed (Njoku et al ., 1980) to account for moisture transport in the soil 
so that it also predicts the spatial and temporal variations of volumetric 
soil water content. This involves a numerical process more complicated than 
the dry model calculation, since it requires the simultaneous solution of a 
pair of coupled, parabolic differential equations for heat and moisture flovf. 

All the aircraft and HCMM data sets have been evaluated to ensure a 
sufficiently good ensemble of data in order to produce a set of high quality 
thermal inertia images suitable for analysis. This includes eliminating cases 
with excessive noise in the thermal or albedo data. Hie daytime thermal, 
nighttime thermal, and albedo data sets were coregistered for each individual 
data set before they were used in the model: once without accounting, for 

topographic corrections, and in the case of the HCMM data, again with the 
corrections included. Topographic Information (for slope azimuth and orient- 
ation), in the form of digitized values from the National Cartographic 
Information Center (NGIC), is added to the proceaslng in the latter case. 
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A . Aircraft Data Processing 


1. M^S Aircraft Scanner Data Processing 

Included in this study were ten sets of flight lines using the 
Bendlx Modu3,ar Multispectral Scanner (M^S). The M^S was mounted in either a 
NP-3A or a l<fC-1.30B aircraft (Table 1). The M^S is an eleven channel scanner, 
with ton narrow bandwidth channels in the visible to neat Infrared range and 
one broad bandwidth thermal channel (Table 5), When operating at night, only 
the eleventh (thermal) channel is activated. 


TABLE 5 

M^S WAVELENGTH BANDS 
Channel Bandwidth 


Number 

(Mi 

crometers) 

Detector 

1 

.33 

- 

.44 

Silicon 

2 

.44 


.48 

Silicon 

3 

.49 

- 

.54 

Silicon 

4 

.54 

- 

.58 

Silicon 

5 

.58 

- 

.62 

Silicon 

6 

.62 

- 

.66 

Silicon 

7 

.66 

- 

.70 

Silicon 

8 

.70 

- 

.74 

Silicon 

9 

.76 

- 

.86 

Silicon 

10 

.97 

- 

1.065 

Silicon 

11 

7.95 

- 

13.5 

HgCdTe 


Fot the thermal channel, there Is an internal calibration to 
teierence blackbodies* Once during each scan line the thermal sensor monitors 
two blackbodies, one warmer, the other cooler than the anticipated ground 
temperature. Tlie actual temperatures of these two blackbodies are indepen'- 
dently monitored by thermistors and automatically recorded with the scanner 
data. 

The channel U DN value for a pixel is assumed to be 

linear with thermal radiance (r^) rather than temperature, but since radiance 
of a blackbody is given by 

tj. “ eor^ (2) 

the temperature can be found from 

A 

- aDN, , f b (3) 

eo ii 

where a and b are determined by solving EqUvatlon 3 for the two blackbodies of 
known T and DNj^j. The emissivity (e) is assumed constant. 

There is no corresponding internal calibration for channels one 
through ten. As a result, calibration of these channels required on-site 
ground measurements, taken at the time of the flight. The calibration used 
for the March, 1975 flights has been previously reported by Gillespie and 
Kahle (1977). The calibration process was repeated for the October, 1976 
flights and is described below. This second calibration was used for all 
subsequent flights. 

Spectral reflectance was measured at seven sites within the 
Cuprite district during October, 1976, using the Jet Propulsion Laboratory 
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Fortablo Field Reflectance Spectronieteir (PFRS) (Goetz et al.^ 1975). Tliese 
seven sites were chosen for their diverse spectral characteristics and for 
their ability to be located on the digital images. 

It was assumed that the response of each sensor is linear with 
respect to the reflectivity of the observed region. Hiat is, for each 
channel n, 


DNn » m„Rn f b„ 


(4) 


where R^ is the average reflectance within the bandpass of channel n. Plots 
of DNj^ vs. for the seven calibration sites are shown in Figure 13. 

Channel I, whose bandpass is outside the range of the PFRS and represents a 
relatively minor portion of solar power, was not calibrated. The coefficients 
m^ and were computed from these data, using a least*-8quares fitting 
algorithm. 

Since the value of interest is the amount of energy absorbed by 
the surface, the desired overall albedo is an average of the channels, weight- 
ed by the amount of solar energy Incident upon the surface within each band- 
pass. That is, the weighted albedo, A^, is represented by 


10 


I 

n*2 


W R 
n n 


w 


10 

1 w 

n“2 


(5) 


where the weighting factor is computed as 

00 

- / F (A)S (X)dX 
n Q n, IT 


( 6 ) 
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In Bkjuatlon 6 is the product of the filter transmittance and scanner 

detector sensitivity (overall sensor response), S^(X) Is the solar irradiance, 
and X is the wavelength. Combining Equations 4 and 5 yields 


where. 


and 


A 


w 


10 

^ I 

n=2 


U DM 
n n 


C 


10 

I 

n-2 


W b 
n n 

m 


10 

I 

n«2 


W 


u 

n 


W 


m 


n 


10 

I 

n=2 


W 


(7) 


( 8 ) 


(9) 


The computed values for C and are listed in Table 6. 

After a day- nigh*; pair of images has been calibrated, the 
images must be geometrical ..,, ectified. Each image was resampled (using a 
linear interpolation) to remove the distortions originating from the panorama 
effect and from a non-urtity aspect ratio. The roll, pitch, and yaw of the 
aircraft carrying the scanner introduce additional geometric distortions. 
While these distortions were sometimes quite noticeable, no attempt was made 
to remove them. 
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In oiHietr to Renerntt* « nlght“tO”(l«y tomp«.'!ratniHi Ulf 


t\w nli’.httirac' ImiiKo munt ho coglstcted to tho dnytlinc tiwigo. J)ay-^ivlght 
roRl.Htration ha» BOim>timos proven quite (.Ul’Meult, ‘n\e prUnary problem In 
ntrcrnft imftge regiatcation arisen from the roll, pitch, and yaw dlntortlouH 
within the Imagea. tliese distortions are often local In nature and completely 


■CAHI-K () 


CAhllVRM'ION 


COEPPIUIKNTS FOR 


SCANNER 


2 

3 

4 

5 

6 


7 

8 
9 


10 


c 


Ji/aiuj>^ 

.000273 

.000333 

.000324 

.000363 

,000324 

,000341 

.000347 

.000502 

.000388 

-.0946 


uncocrelated from day to night. A secondary registration problem Is that the 
features visible In a nighttime thermtul Ittuige are not necessarily conspicuous 
In a daytime Image, visible or thermal. The techniques used to determine a 
registration transformiitlon fall into two general classes: 1) polynomial 

fitting, and 2) trlangvilatlon. 


Xn a poXynomlnl fit, the ttansforiwition from the old coordin- 
ates to new coordinates Is described by a simple polynomial in two variables 
(line and sample). Witli a first order polynomial fit, even a slight deviation 
from a straight flight line yields substantial misregistration of most of the 
itTOge. Higiicr order polynomial fits (up to third order have been tried) are 
somewluvt better at tracking non-linear and local distortions, but tend to work 
very poorly near the edges and corners of an image. 

In a triangulation scheme, the image is divided into a grid of 
numerous rectangular cells. The transformtition for each cell is a linear 
transformation, but the transformation is independently recomputed for each 
cell, using only the nearest tiepoints that enclose the corners of the cell. 
This technique yields a generally closer registration than the polynomial fit, 
hut there imiy be unpleasant artifactual features along cell boundaries. Tl^ere 
is no guarantee that the transformations are continuous across cell boundaries 
and, as a result, such a boundary may appear as an apparent tear, crimp, or 
fold in the transformed image. Tlie triangulation method may also lead to poor 
registation of corners and edges (though not to the extent of the high-order 
polynomial fits). Most images were registered using the triangulation 
technique. 

2. HCM and DMS Aircraft Scanner Data Processing 

Computer compatible tapes of nine sets of day-night flights 
using scanners mounted on a U-2 aircraft were received (Table 1). Seven sets 
utilized the HCM scanner, while two sets used a modified Daedalus model DEI- 
1260 multispectral scanner (DMS). 

The HCM scanner simulates the two sensors used in the HCMM 
satellite. Tl\e DMS scanner contains eleven channels, ten with visible and 
near IR bandpasses, and one channel in the thermal range. The characteristics 


of these two scanners are 8unui»r:l/,e<l in Table. 7. Xn all cases, the calibra-' 
tion Infornwitlon was pvovlded by the operators of the scanner. A more detail- 
ed description of the HCM scanner characteristics may be found in Murphy _et^ 
al. (1978). 


TABLE 7 

CHARACTERISTICS OF U-2 BORNE SCANNERS 





ncM 

DMS 

Nominal 

Altitude 


65,000 ft. 

65,000 ft. 

Angular 

Resolution 

(lEOV) 

2. 8 mr 

2.5 mr 

Spatial 

Resolution 

(at Nadir) 

55 meters 

50 meters 

FOV 



90® 

85® 

Samples 

Per Scan Line 

400 

715 

Biindpassea 

1 

.10,5-12.5 Hiv 

0.38-0.42 



2 

0.51-0.89 

0.42-0.45 



3 


0.45-0. 50 



4 


0.50-0.55 



5 


0.55-0.60 



6 


0.60-0.65 



7 


0.65-0.69 



8 


0.70-0.79 



9 


0.80-0.89 



10 


0.90-1,10 



11 


10,40-12.50 
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Noise (both coherent and random), missing segments of data, 
severe geometric distortions, and cloud cover combined to make most of these 
day-night pairs unusable (Figure 14). Only the March, 1978 and April, 1979 
data of Death Valley were adequate for generation of thermal inertia images. 

The 1:250,000 NCIC digital elevation data were used as the 
basis for registering daytime and nighttime Images (Figure 15). With the 
images associated with elevatiou data, topographic data could be used to gain 
a more sophisticated estimate of the thermal Inertia. (At a resolution of 
63.4 meters, the NCIC data are comparable to the U-2 images, but Inadequate 
for the higher resolution images from the M^S scanner). 

Tlie registrations were performed in the following sequence: l) 

the nighttime image was registered to the daytime image; 2) the daytime albedo 
inwage was used to find a transformation to the elevation image; 3) this 
transformation was applied to both the day and night (registered to day) 
Images. 

This two-stage registration of night to map base Implies a 
summation of errors from each transformation. The alternative, a direct 
registration of the nighttime image to the NCIC data was less satisfactory, 
due to a lack of distinct features common to both images (compare Figures 15b 
and l6c). 

B . HCMM Satel llte Data Processing 

Nine sets of HCMM overpasses were Included in this study 
(Table 2). In general, a set consists of two nighttime scenes, one day apart, 
and two daytime scenes, also one day apart, with the first daytime scene 
following the second nighttime scene by twelve hours. That is, the scenes are 
in the order night-night-day-day (NNDD). In several cases, the test site data 
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Figure 14a. Fxamiile of cohereiiL iiol.se ami .si'vere geometric distortion. 
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Figure 14b. Example of random noist 
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Figure 14c. Example of missing data segments. 
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Figure 1.5a, Shaded topographic relief image of Walker Lane, Nevada. 
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Klgurt' 151). Shack'd t opo>;rai)hic reltc'f iiniigc* of Doafh Valloy, California. 
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Figure 15c. Shaded topographic relief image of Pisgah, California. 
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Fijiuri.' 16a. Death Valley March 1978 registered visible image 
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Figure Ihc. Death Valley March 1978 registered night thermal image. 
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were contained on two tapes; these data were mosaicked together before any 
other processing was performed. 

The characteristics of the HCMM satellite and its sensors are 
reported in detail elsewhere (Bohse, et al . , 1979). The data received a.t JPli 
were calibrated data for both channels, and geometrically rectified to a 
Hotine Oblique Mercator base. The spatial resolution was approximately 500 
meters. Registration was a multi-step operation, and involved the following 
steps for each set: 

1. For the entire scene, both nighttime images and the second 
daytime image were registered to the first daytime image of the set. In a 
typical registration, about 25 points common to both triages would be located 
and used to find a geometric tranformatlon. The transformation, utssd was a 
first order, unconstrained polynomial, i.e. , 


x' = Ax f By f C 


y' = Dx 4- Ey f F 


( 10 ) 


where the coefficients A, B, G, D, E and F are determined by a least-squares 
fit. The average residual to the fitted surface was usually under 2.0 pixels. 

2, For each of the test sites, a day-night pair was chosen for 
generation of thermal inertia Images. Whenever possible, a pair separated by 
12 hours was chosen, but if a site was missed by one of those passes, a 36 
hour pair was chosen. 

3. The local area surrounding each test site was re-registered, 
night to day. These corrections were limited to simple offsets, i.e.. 


/ 
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X' - X f A 
Y' “ y f B 


(U) 


The relative scale and rotation of the images were not allowed to change for 
these smaller scenes, since the overall registration was expected to yield a 
better estimate of these values. Offsets varied from 0 to 2 pixels. 

4. Except for the. San Rafael pair, for which topographic 
information was not readily available, all pairs were then registered to the 
1:250,000 NCIC digital elevation data. A first order unconstrained fit was 
used, with typical residuals of about one pixel. 

As a result of this technique, all pairs were registered to a 
common Universal Transverse Mercator (UTM) map base at each site. 

C. Generation of Tl'iermal Inertia Images 

A model has been developed at JRL for the purpose of predicting the 
temperature of the Earth's surface as a function of albedo, topographic 
orientation, and meteorological conditions. Using this model at each site and 
date, day-night temperature differences were predicted for a range of albedos 
and (when topographic information was available) ground slopes and slope 
azimuths. These predicted AT's were then used to generate a table that ex- 
pressed thermal inertia as a function of albedo, AT, slope, and slope 
azimuth. To generate a thermal inertia image, each pixel was referred to this 
table and assigned the value that corresponded to each of those variables. 

The theoretical basis and method of application of this model have 
been described in detail elsewhere (Kahle, 1977). In brief, the model works 
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In the following manner. An explicit finite difference scheme is used to 
approximate a solution to the one-dimensional heat flow equation 


at' 



( 12 ) 


where 'T is temperature, k is thermal dlffuslvity, t is time and z is depth. A 
one centimeter vertical grid is used from the surface to a depth of z = 50 
centimeters, and temperature values are recomputed at each depth for each 
twenty second time interval. Initial temperature profile values are estimates 
(derived from on-site measurements when available); the importance of these 
values is minimized by initializing the model to several hours prior to the 
time interval of Interest. The lower boundary condition is that the 
temperature at 50 cm depth is constant. The upper surface boundary condition 
is that the heat flux is continuous across the boundary. 

In order to compute the fluxes at the surface, some meteorological 
and geographic data are needed. The date, latitude, and surface elevation are 
needed to com|>ute solar radiation. Air and ground temperatures and wind speed 
are used to compute the sensible heat flux. If the latent heat flux is 
Included, the. air and ground specific Jiumidltles are also used. 

Thermal Inertia images were then created for each date at each site, 
with and without the topographic corrections. Some selected full frame 
thermal Inertia Images were also produced. Histograms of the distribution of 
numerical values of thermal inertia in each scene were also generated during 
this process. Examples are discussed in subseauent sections. 
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V. DATA AND TECHNIQUE EVALUATION 


We have worked entirely with digital HC^tM data in an effort to produce 
quantitatively meaningful thermal inertia images. In this section we discuss 
and assess many of the elements involved in a numerical study of the data, the 
problems encountered, and the resulting thermal inertia products. Specif ical“ 
iy we discuss atmospheric transmission and emission effects, problems of HGMR 
data calibration, the dependence of thernuil inertia calculations on meteor- 
ological and soil moisture conditions, and model sensitivity studies. We 
discuss the numerical values of thermal inertia which we obtain, as a function 
of site, material and time, and finally, we compare modeled thermal inertia 
With apparent thernuil inertia. 

A. Atmospheric Effects 

In order to determine the temperature of the surface of the Earth 
from radiance measurements obtained with a satellite sensor, there are several 
factors which must be taken into account. These Include the emlssivity of the 
surface, atmospheric attenuation of the upward radiation, atmospheric emission 
of radiation upward into the sensor and downward to be reflected upward to the 
sensor, and finally the calibration of the sensor, relating received radiance 
to apparent brightness temperature. Hie sensor calibration is discussed in a 
later section. We will discuss the. other factors here. 

The spectral radiance received by the satellite sensor can be ex- 
pressed by 

Lg(X) - [E;^LBga,Tg) f (l-e^)L^^(X))T;^ (13) 

Here the ground temperature, T„, is the quantity we wish to determine, from 

O 

the measurement of Lg(X). Therefore it is necessary to measure, model, or 
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Ignore the rest of the terms in Equation 13, as defined and discussed below. 

llie spectral radiance emitted by the ground is given by 

Lg(X) “ Lgg(^,Tg) (lA) 

where is the emlsslvity of the surface material and liBg( black” 

body spectral radiance at wavelength X and temperature T. In the 10.5 to 
12.5 viti range, most natural surface materials have emissivlties between about 
0.90 and 0.98. Our lack of knowledge of the emlsslvity of the surface 
material will thus Introduce an error of a few percent if we assume an average 
value around 0.95. 

The reflected portion of the downward atmospheric radiation 
L^|(X) is small because the reflectivity of the surface (l~Ej^) is small, on 
the order of 0.05, and the atmospheric downward radiation L^^(X) is also small 
in this window region compared to the radiance from the surface. Therefore, 

.iii can ignore this term. 

The remaining quantities which we need to determine from models, 
measurements or some combination of these, are the atmospheric transmission, 
and the upward atmospheric radiation L^^(X). These two terms are related. 
In these wavelengths, contributions to both are due almost entirely to the 
presence of water vapor and aerosols. As the concentrations of these absorb” 
ers-emitters and scatterers increase, the contribution of the ground signal to 
the radiation received at the satellite will decrease and the contribution of 
the atmospheric signal will increase. We turn to the theory of radiative 
transfer in order to understand the effects. 

The basic problem of radiative transfer is to determine how a radia- 
tion field is altered on passing through and interacting with a medium 
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(Ohttndrasckhor, 1960). If wo considor a poacll of radiation of intensity ly 
(in a frequency Interval V to v I dv) pcopaj^ntlng tiirough a roedluia In a sped'* 
fled direction, then we can consider the various ways this Intensity will be 
altered In travellnj^ a distance ds. 

Tlie radiation c.an be weakened by two processes: true absorption by 

the medium; and, scattering of the radiation Into another direction. In 
practice, both of these processes arc combined into a mss attenuation coeffi- 
cient, defined as Ky In 


dly — (CyplydS 


(15) 


where p Is tV'.e mass density of the inedlura. This attenuation coefficient can 
be due to true absorption only, scattering only, or a combination of both. 

As with the weakening of our pencil of radiation In traversing an 
element of mass, there are two similar processes enhancing the radiation: 
scattering of radiation into the direction of the pencil, and emission by the 
mass element. Again, the two are treated together, this time with a single 
emission coefficient jy defined such that the element of mass dm emits Into 
the solid angle dw an amount of radiation In the frequency range v to v h dv 
equal to jy dm dw dV in unit time. 

Tlie source function is defined as 


J 


V 



(16) 


For scattering only, this is given by 
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4tr 




w 2 tt 

/ / )stn 0^ d0' df 

0 0 


U7) 


where p Is the phase function of scatterlrg, when we are considering the scat- 
tering from direction 0', 4>" , into direction 0, <(>. 

For no scattering and an atmosphere In thermal equilibrium, the 
emission coefficient according to Klrchhoff's Law can be given In terms of the 
absorption coefficient 


where By(T) Is the black body radiation. Thus the source function J ^ for 
atmospheric emission Is just B^(T). In general the source term will be the 
sum of the scattering and emission source terms. 

The basic equation of radiative transfer sums up the contribution of 
the various processes which we have just discussed 





(19) 


or 


dl, 

ds 


tc pi - tc pJ , 


( 20 ) 


The solution to the equation of radiative transfer is 


I(s) - I(o)e"'^^®»°^ f / J(s')e“'^^®’®'^<pd8' 


( 21 ) 
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where now t is the optical thtckncss defined by 


8 

t(s,s') “ / Kpd8 (22) 

s' 

Thus the intensity of radiation at the satellite, I(s), is given by the inten- 
sity at ground level, I(o) , reduced by plus the radiation from the 

atroospheric sources reduced by While this solution is easy to 

write, the evaluation of the atmospheric absorption, scattering and emission 
parameters is not easy. 

In the middle infrared wavelength region in which HCMR operates 
(10.5-12.5 pm) both scattering and absorption-emission occur. Rayleigh scat- 
tering occurs when the particle size is small compared to the Wavelength. 
Hence, Rayleigh scattering from the atmospheric gas molecules is primarily 
confined to the visible pt>r'tiow of the spectrum. Hie scattering in the middle 
Infrared wavelengths of concern here Is by aerosols — liquid water, ice and 
other particulates. 

The major atmospheric constituents, N 2 and O 2 , being symmetrical 
homonuclear diatomic molecules, without a permanent dipole field, do not have 
vibrational and rotational absorption features, and are essentially trans- 
parent to visible and infrared radiation. The absorption and emission in the 
middle Infrared wavelengths is caused by the minor polyatomic gases In the 
atmosphere, primarily H 2 O, CO 2 and 0-j. Even though the 10.5 to 12.5 pn region 
is considered a g id atmospheric window, there is significant absorption and 
emission by the water vapor. Unfortunately, the water content of the atmos- 
phere is quite variable and the radiation effects are extremely sensitive to 
these changes. Contrary to the statements in the HCMM Data Users H andbook 
(Goddard Space Flight Center, 1978), these effects are much more sensitive to 
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water vapor amount than to atmospheric temperature, as can be inferred from 
the source term in the radiative transfer equation. 

Various numerical solutions to the problem of atmospheric radiative 
transfer exist. Ilhey vary in complexity and accuracy from those giving a 

single value of total radiation based on one or two average atmospheric param^ 

eters, to numerical models which follow every known absorption line, as a 
function of measured pressure, temperature and atmospheric gas profiles. For 
the HCMM problem we required a model of intermediate complexity. 

A good example of a calculation of the radiation emerging from the 
top of the atmosphere, developed to study the atmospheric effects on the 
satellite measurement of surface temperature can be found in the work of 
Anding e.t al. (1971). One of their figures, r^^produced here as Figure 17, 
illustrates many of the important effects. The atmospheric model shown in 
this figure represents a mean atmosphere for 30% latitude. The solid curves 
show the spectral radiance emergir- , -‘nm the top of the atmosphere as a 

function of surface temperature. The dashed curves show the black body 

radiation which would reach the satellite if there were no intervening 
atmosphere. The solid curve labelled 0% represents the atmospheric component 
of the upward radiation. The magnitude of atmospheric effects in the 10.5 to 
12.5 un wavelength region of interest is readily apparent. Because this model 
was developed for studies of sea surface temperatures, it assumes the surface 
to be at sea level. 

Another model, which we have been using in our numerical studies. of 
HCMM data and which gives results similar to those of Anding et al. (1971), is 
the LOWTRAN 5 model, developed at the Air Force Geophysical Research Labor- 
atory. The source code for this versatile computer model is available to the 
public. The model includes a choice of six standard atmospheres: U. S. 
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Figure 17. Upward spectral radiance at 100 km, with surface temperature as a 
parameter, after And in g ct (1971). 
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Standard Atmosphere, tropical, mldlatttude summer, midlatitude wlntetj sub- 
arctic summer and subarctic wlnter> Alternatively, it also allows the user to 
enter his own atmospheric model or measurements* Figure 18, after Knelzys 
al* (1980), shows some LOWTRAN 5 results. The top figure Illustrates the 
extreme sensitivity of atmospheric transmission to the particular atmospheric 
model selected. Note that the atmospheric transmittance in the HCMR wave- 
lengths varies from approximately 50% to 95%, pritmitlly as a function of water 
vapor content. Another important feature of the LOWTRAN 5 model Is that it 
allows the surface elevation to vary. In our running of the model over the 
range of elevations present at our test sites (- 73 m at Badwater in Death 
Valley to 3400 m at Telescope Peak only a few ral.les away). It became apparent 
that meaningful atmospheric corrections must include surface elevation. 

In order to determine the magnitude and variability of the atmos- 
pheric corrections required for HCMM data interpretation, we have made several 
runs with LOWTRAN 5. Assuming a given surface temperature, elevation and 
atmospheric model, we used LOWTRAN to compute the atmospheric transmission, 
and the upward radiance at the top of the atmosphere from both the ground and 
atmospheric emission. This radiance was converted to the effective temper- 
ature which would be sensed by the satellite, and the difference between 
surface temperature and effective temperature was noted. Also the peccentage 
of upward radiance originating in the atmosphere rather than at the surface 
was calculated. We then allowed the parameters to vary, one at a time. Ttiese 
results are presented in Tables 8-12. 

In Table 8 we used the rawinsonde data for July 18, 1978 at 2300 Z 
(1500 PST) and an assumed ground temperature of 300^ and allowed ground 
elevation to vary from X km to 3 km. The atmospheric effect was to make the 
ground appear colder than it actually was by from 4°K at an elevation of 1 km 
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Figure, 18. (a) transmittance, and (b) 

radiaiice spectra for a VGrtical 
path looking at the ground from 
space, for the six model atmo- 
spheres of LOWXRAN 5, aft.er 
Kneixys (1980) . 
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to only at an elevation of 3 km. In Table 9 we show the result when 
elevation is constant at I km, and the surface temperature is allowed to vary 
from 260^K to 340“IC. Tliere is a striking variation in the atmospheric effect, 
with the apparent temperature of the surface varying from 3‘^K warmer for the 
coldest surface to ll°K colder than the actual surface temperature for the 
hottest surface, 'finis strong surface-temperature dependence of the atmospher- 
ic effects is in good agreement with effects noted by both European investi- 
gators and by us in comparisons of measured surface temperature and satellite 
data. Table 10 shows the results for several different rawinsonde 
observations over an eight day period In July, 1978; ground elevation and 
temperature are held constant. Only very small differences are noted. 

Finally in Tables 11 and 12, we repeat the conditions of Tables 8 and 9, only 
now using the LOWTRAN midlatitude summer atmospheric model. This model, 
unlike the rawinsonde data, is available down to sea level. Tliiere is somewhat 
more atmospheric water vapor in this model than in the rawinsonde observations 
So the effects are slightly larger, but quite similar to those noted in Tables 
8 and 9. 
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'CABLE 8 

results from lowtran using rawinsonde data for varying surface elevations 


Atmospberjlc Model 

a.irfrtce 

Blevatlon 

S^irfacp 

Tempernturs 

AtmoBpherie 
Tiranemj .flalon 

% Rullanet*. 
from Armofi. 

Effect;! VP, 
Temperature 

5T 

RAwl'isonde 

July 18* 1978 1,500 PST 

LO 

300 

.7468 

20,00 

296 


M )l 

1.5 

It 

,8l3f 

13.98 

297 


H M 

2.0 

II 

.8650 

9.53 

297 

-3 

II It 

2.5 

II 

.9015 

6.61 

298 

-2 

II It 

3,0 

It 

.9272 

4 , 59 

299 

»1 


TABLE 9 

RESULTS FROM LOWTRAN USING RAWINSONDE DATA FOR VARYING SURFACE TEMPERATURES 


Atraoepheric Model 

Surface 
R1 ovation 

Surface 
Tempo rat are 

Atmospheric 

Transmission 

% Radiance 
from Atmos. 

Effective 

Temperature 

61.’ 

Rfiwlnsonde 
July 18, 1978 

1500 PST 

,l 

260 

.7468 

31.19 

263 

3 

II 

tl 

It 

270 

II 

27,66 

271 

1 

H 

It 

H 

280 

II 

24.72 

280 

0 

It 

It 

II 

290 

tl 

22.14 

288 

-2 

M 

II 

II 

300 

11 

20.01 

296 

-4 

It 

It 

0 

310 

II 

18.16 

304 

-6 

It 

II 

II 

320 

tt 

16.60 

312 

-8 

II 

II 


330 

It 

15.20 

320 

-10 

It 

If 

tl 

340 

It 

14.07 

329 

-11 
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TABLE 10 



RESULTS 

FROM LOWTRAN 

USING RAWINSONDE DATA 

FOR VARIOUS 

DATES 


Atmospheric 

Model 

Surface 

BlevaCton 

Surface 

TeKperature 

Atmospheric 

Transmission 

% Radiance 
from Atmos. 

Effective 

Tempuroturfi 

6T 

Rawinsonde 
July 16, 1978 

0300 

PST 

1 

300 

.7463 

20.04 

295 

-4 

July 

16, 1978 

1500 

PST 

It 

II 

.7463 

20.03 

It 

II 

July 

17, 19/8 

0300 

PST 

M 

II 

.7453 

20.14 

tl 

II 

July 

17, 1978 

1500 

PST 

M 


.7471 

19.97 

II 

If 

July 

18, 1978 

0300 

PST 

ft 

l( 

.7471 

19.98 

II 

II 

July 

18, 1978 

1500 

PST 

M 

II 

.7468 

20.00 

It 

M 

July 

21, 1978 

030-0 

PST 

ft 


.7485 

19.89 

11 

II 

July 

21, 1978 

1500 

PST 


It 

.7459 

20.11 

It 

II 

July 

22, 1978 

0300 

PST 

'll 


.7478 

19.96 

II 

If 

July 

22, 1978 

1500 

PST 

tl 

It 

.7469 

20.05 

II 

It 

July 

23, 1978 

0300 

PST 

il 

II 

.7476 

19.93 

II 

II 
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TABLE 11 

RESULTS FROM LOWTRAN USING MIDLATITUDE SUMMER MODEL FOR VARYING SURFACE ELEVATIONS 


Atmospheric Model, 

Surfocc 
Ele.va tlou 

Surface 

Temperature 

Atmospheric 

Transmission 

% Radiance 
from Atmos. 

Effective. 

Temperature 


Mid loti tilde Summer 

d.o 

300 

.5627 

38,37 

295 

-5 

M M 

0.25 

It 

.6167 

32.93 

295 

-5 

II II 

0,50 

It 

.6657 

28.07 

296 

“4 

M M 

0.75 

II 

,7096 

23.76 

296 

-4 

II It 

1.0 

II 

.7488 

19.95 

296 

-4 


1.5 

It 

.8144 

13.97 

297 

-3 

M II 

2.0 

II 

.8638 

9.55 

298 

-2 


2.5 

II 

.8986 

6.70 

298 

-2 


3.0 

M 

.9230 

4.72 

299 

-1 


RESULTS FROM LOWTRAN 

TABLE 12 

USING MIDLATITUDE SUMMER MODEL FOR 

VARYING SURFACE TEMPERATURES 

Atmospheric Model 

Surface 
Eievat ion 

Surface 

Temperature 

Atmospheric 

Transmission 

X Radiance 
from Atmos. 

Effective 

Temperature 

SI 

Midlatitude Summer 

1 

260 

.7488 

31.11 

264 

4 

II If 

H 

270 

II 

27.57 

271 

1 

II II 

II 

280 

II 

24.63 

279 

-1 

II 11 

It 

290 

II 

22.09 

288 

-2 



300 

II 

19.95 

296 

-4 

II It 

II 

310 

It 

18.10 

304 

-6 

II II 

II 

320 

It 

16.54 

313 

-7 

II *1 

II 

330 

II 

15.20 

320 

-10 


II 

340 

II 

14.04 

327 

-1 1 


II 

350 

H 

13.01 

338 

-12 
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Elevation data are available for the test sites In digital format, 
and all the rawinsonde data are on hand. However we have not had time to 
develop the computer algorithms to apply the atmospheric corrections differ- 
entially across the test sites. This has been left for a follow-on study. 

B. HGMR Calibration 

Our processing of the HCMM thermal data Included an experiment 
designed to help verify the satellite calibration. Tlie effective temperature 
corresponding to the radiance at the top of the atmosphere is given by the 
formula for infrared data conversion supplied by the HCMM Project Office at 
GSFC, 


T^(DN) 


K2 


In [ 


K1 


DN-K3 


■f 1 


(23) 


in which DN is the digital number value (0-255) of a pixel; T„ is the temper- 

ature in degrees Kelvin, with T(0) = 260. 0°K and T(255) = 340.0*^K; and 

K1 = 14421.587, K2 = 1251.1591, K3 = -118.21378. T, differs from the surface 

s 

temperature T due to the atmospheric effects discussed earlier. Post-launch 

o 

HCMR calibrations were performed by the HCMM Project Office, based on compar- 
isons of the satellite-measured radiance (and hence temperature) and the 
"effective temperature" at the top of the atmosphere calculated from ground 
measurements at White Sands, New Mexico and a GSFC atmospheric correction 
program. Based on such calculations during the first sixty days after HCMM 
launch, a f5.2°K offset was noted, with the HCMR data warmer than the surface 
data. The project office decided to adjust the calibration constants to 
offset all data by f5.2°K, to match the surface "truth," and investigators 
were so notified. A second "validation" calibration was undertaken by GSFC a 
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few months later, using the same technique, and again using White Sands data, 
and with the HCMM data now offset by 4-5. 2®K. This time a 5.5®K. offset in the 
other direction was noted (or actually only a -0.3^K offset from the original 
calibrations). However, the offset was now being routinely added to 

all data. At the time that investigators were informed of this problem by 
GSFC personnel, they also discussed two other possible sources of error due to 
noise and drift of up to 0. 4°K and l.O^K respectively. Finally the HCMM 
project office sent a memo stating that both of the above calibration studies 
were based on their incorrect atmospheric model. Their latest estimate, using 
their new atmospheric model, showed that the above offsets should have been 
4-4. 0°K and -5.7‘'’K (instead of f5.2'^K and -5.5°K). Meanwhile, several European 
investigators were stating that they found surface temperature dependent 
offsets in the data tanging from about 5 to 10°K, increasing with higher 
surface temperature. 

The above if formation continued to arrive from GSFC during our 
attempts to put our data analysis on a sound quantitative basis. Recognizing 
that tl -are appeared to be several uncertainties in the data, we first chose to 
compare effective temperature measured by the satellite with our measured 
surface temperatures, with no atmospheric corrections. From this study we 
hoped to establish at least an estimate of the magnitude of the calibration- 
offset and the atmospheric correction problems. Our goals in this experiment 
were to estimate the size of this temperatvire difference and to determine 
roughly how much it varied from day to day. 

For this investigation we used ground truth data from two sources. 
The primary source was field data taken by JPL personnel near Goldfield, 

Nevada in July, 1978. Temperatures were meiffsured with the Precision Radiation 
Tliermometer, model PRT-5 of the Barnes Engineering Company, at eight different 
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sites at approximately the time of the satellite overpass. These sites were 
located on the HCMM albedo image with the aid of air photos, and DN values 
were taken from both day and night thermal images which had been registered to 
the albedo image. Temperatures computed from these DN values using the con- 
version formula, Equation 23, were then compared with the recorded temper- 
atures . 

In order to compare measurements over a wider temperature range, 
recorded temperatures were obtained from additional sources. National Weather 
Service offices in Las Vegas, Nevada and Los Angeles, California provided 
shoreline water temperatures at, respectively, Lake Mead and six beaches in 
the Los Angeles-San Diego area. Offshore water surface temperatures for Lake 
Tahoe were obtained from Professor Charles R. Goldman of the University of 
California at Davis. Ground surface temperatures were provided by the Nation- 
al Park Service at Death Valley National Monument in California and by the 
United States Army at the Yuma Proving Ground in Arizona. Measurements from 
the latter two sources were subsequently excluded from the test as, on the 
dates under consideration, the temperatures recorded were outside the range of 
the HCMR. An attempt was made to include the melting temperature of snow as 
additional ground truth data but the resolution of the HCMM images was too 
poor to locate accurately the edge of any snow masses. 

We sought a functional relationship between the satellite measure- 
ments, Tg, and the temperatures recorded at the surface, Tg, using a linear 
least squares analysis. Several data pairs, with x = Tg and y = Tg were used 
to compute the constants a and b which would give the best least squares fit 
to three different functions: linear, y = a f bx; power, y = ax^; and expo- 
nential, y = ae^^ (or logarithmic, y = a f blnx). The result, Tg = Tg(Tg), 
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which was obtained from the best fit tc the data was then compared with the 
assumption Tg = Tg, which pertains to a dry atmosphere. 

Initially T„-T„ pairs for three different dates were examined 
separately as shown in Figure 19. The difference in T„ “ T„(T„) between the 

b b g 

two dates, July 17, 1978 and July 22, 1978 - both clear, windy days as 
observed during the acquisition of the field data near Goldfield - ranged from 
0 to 3°C. For both dates, the satellite-derived temperature was lower than 
the recorded ground temperature. No curve is plotted for the 19 August 1978 
data because the fit to each of the three trial functions resulted in a very 
low correlation. These data points are in good agreement with the other two 
curves. Because field data, the primary source of data used in this calibra- 
tion experiment, were available only for summer months, possible seasonal 
effects were not examined. 

Because the daily variation in T^ was small we combined all the data 
pairs for temperatures measured during the day. Tlie result is shown in Figure 
20. A corresponding result for the nighttime data pairs could not be derived 
due to poor correlation of the data. A final comparison of Tg and Tg was made 
by combining all the data pairs. As shown by the dashed line in Figure 21, 
the addition of the nighttime temperatures did not greatly affect the form of 
'^o “ To(To)‘ is seen to be 5-6°C lower than T near 20®C and 8-9®G lower 
at higher temperatures of 50-60°C. These differences are fairly consistent 
with the atmospheric effects we calculated in the previous section. The 
offset for lower temperatures agrees with that found by the HCMM Project 
Office in their second calibration experiment in October 1978 although they 
had Included atmospheric effects in their calculations, and the temperature 
dependence agrees with the effects noted by the European investigators. 
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Figure 19. Best fit curve for Tg-Tg data pairs for different dates. The 
ideal case is shown by the line labeled Tg = Tg. 
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POOR QUALITY 



Figure 20. Best fin curve for daytime Tg-Tg data pairs for all dates. Nighttime 
data pairs are denoted by crosses. The ideal case is shown by the 
line labeled - Ta . 






After receiving notification by OSFG on Augunt 18, 1980 of the error 
in the radiative transfer program used in their calibration we decided to use 
only the formula for Infrared data conversion, Equation 23, in our processing 
for this study. It is our intent in a follow-on study to use our atmospheric 
programs, along with digital elevation information and the ground data dis- 
cussed above, to better test the HCMR calibration. 

C. Dependence of Thermal Inertia Calculations on fletoiorological 

Conditions 

Heat and mass fluxes between the ground and air are required when 
calculating the thermal inertia of surface materials. In particular, sensible 
and latent heat flux densities are often the predominant terms in the equation 
for the surface energy balance. Moreover, evaporation from moist soils 
changes the thermal inertia of the near surface soil layers with depth and 
time. 

The problem of computing sensible and latent heat flux densities for 
use In thermal inertia calculations requires a knowledge of the mesoscale 
variability of meteorological conditions over the area covered by the Image, 
on the order of tens to hundreds of square kilometers. Such detailed inform- 
ation is prohibitively costly to obtain and process, even for a single image, 
since it would require an observational network on a scale of about one 
station per square kilometer. Furthermore, theoretical analysis of the 
variation of turbulent heat flux densities over heterogeneous surfaces is 
still relatively undeveloped. The usual and only practical way around these 
difficulties is to assume that one or two point measurements are sufficiently 
representative of the area. This is a fairly good assumption when surface 
winds are moderate to strong, i.e. more than about 5 ra/s, because stro'Rg 
mixing of the air tends to produce uniform conditions everywhere. On the 
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other hand, large differences in svirface roeteorological fields can exist on a 
microscale when surface winds are light, i.e., less than about 1 m/s. Despite 
this complexity, sensible and latent heat flux densities are often calculated 
from parameterization methods that employ observations obtained at one point 
in an area (cf. Bhumralkar, 1976). In the following, we describe the methods 
used to calculate the turbulent heat flux densities at the ground for use in 
the JPli thermal model. We Include discussions on the details of the formu- 
lation, the types of meteorological data used in tl.e calculations, and the 
limitations of the methods. 

In this study, the turbulent heat flux densities are calculated from 
either of two methods; 1) the bulk aerodynamic method, or 2) the profile 
method. Ihe chief difference between the two methods is that the first 
employs meteorological data obtained at a single height above the ground, 
while the second employs data obtained at two or more heights. 

The bulk aerodynamic method also requires that the ground temper- 
ature and moisture content be known. For the bulk aerodynamic method, the 
sensible heat flux density (H) is given by the expression (Kahle, 1977), 


H = P^CpCoWd^-Tg) 


(24) 


Where is the surface air density, Cp is the specific heat of dry air at 
constant pressure, Cp is an altitude dependent drag coefficient [= 0.002 f 
0.0006(2/5000) ] y Z is elevation in meters, W is a wind speed factor corrected 
for gustiness (wind speed plus 2 m/s), is the air temperature at a height 
just above the ground (~ 2 m), and Tg is the ground temperature. The latent 
heat flux density (L) is given by the expression. 
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The niajor problem in calculating Ii is that q is difficult to determine. (T 

K P 

for the calculation of H can usually be estimated fairly accurately from 
remote sensing data.) One method of computing q is to assuvw that it is a 

O 

fraction (Mg) of the saturation specific humidity “ 0.622 Cg/P), Mg is 

a ground moisture factor that ranges from 0 for completely dry soil to 1 for 
saturated soil. Wl^en moisture transport in the soil is Included in the 
thermal model, Mg is obtained from the model calculations as the ratio, 
®/6saf 5 Is the volumetric soil moisture content (cm of water/cm of total 
soil) at the surface and ®sat is its value at saturation. ®sat has a fixed 
value for each type of soil (Clapp and Hornberger, 1978). If moisture trans~ 
port is not included in the model, 6 may be computed from a rate equation 
along the lines developed by Deard('""ff (1978). 
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The profile method Is n menns of GalculLfifting turbulent heat fluxes 
at the ground from the mean values of simultaneous roeasureroents of wind speed > 
temperature) and humidity made at several levels in the atroospherie surface 
layer. Tlie surf ace layer is the lowest part of the atmospheric boundary layer 
and usually extends to about 10 m (at night) or 100 m (at mld'-af ternoon) above 
the surface* Here the vertical fluxes of momentum, heat, and water vapor are 
assumed to be independent of height so this region is also called the constant 
flux layeri A portable micrometeorologlcai system has been developed at JPL 
(Kahle e_t^. , 1977) to obtain measurements in the constant flux layer over an 
8 meter altitude. Seven levels of instrumentation - measuring wind speed, 
dry~*bulb air temperature, and wet-bulb air temperature - are spaced loga- 
rithmscally (base 2) in the vertical direction from 1/8 m to 8 ro. Tlie geomet- 
ric mean height (GMH) of the array is 1 m. Tlie data are recorded continuously 
and syeraged over 16 minute intervals. 

In the profile method, the similarity hypothesis of >fonin and 
Obukhov (cf. Montn and Yaglom, 1971) is used to calculate the heat flux 
densities. It describes the vertical gradients of wind (U), temperature (T) 
and humidity (Q) in the constant flux layer as functions of a set of non- 
constant, nondimensional parameters and constant scale fa»-'tor8>, 

* (O 

ftT 


i2.= 

9Z 
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vihere 'h is height above the auirfacej k Is von Karman's constant (« 0.35), 
and I'g are the nondlmensional parameters, and U*, T*, Q* are the constant 
scale factors. TChe parameter C is equal to Zih where L is the Monin - Obukhov 
scale length - which is approximately equal to the vertical thickness of the 
constant flU 7 /S layer. Ordinarily, t is < 0 during the day (unstable condi- 
tions) and C > 0 at night (stable conditions). The surface layer Is called 
unstable when the ground is warmer than the air above it so that the temper- 
ature of the air decreases with height in the surface layer, and stable when 
the opposite condition exists. 1 can be obtained from meteorological observa- 
tions because % is proportional to the bulk Richarson number (Rig) for the 
height interval (cf. Yaglom, 1977); viz,, itt our studlcB the height interval 1 
8 m and Rig is computed at the GMH (» 1 m) of the observations. Finally, h*, 
T*, Q* are calculated by the least-squares fit of the observations to the 
Integrated forms of Equations 28. (The parameters ij^, <|^, and (jfi^ are given as 
a function of ? in Table 13.) Then, H and L are calculated from, 


H - PaCpU*T* 

(29) 

h “ Pa • 

(30) 


Figure 22 is an example of H and L calculated from profile data obtained at 
Death Valley in January, 1977. 

How representative are the turbulent heat flux densities calculated 
from thewe two methods? Tlie bulk aerodynamic method Is often used to provide 
surface values for use in atmospheric general circulation models (e.g. Gates 
^£l^. , 1971). The profile method is more often used to determine momentum, 
mass, and heat transfer on level fields (e.g. Morgan et al ., 1971). Ostens- 
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JANUARY 17, 1977 


Figure 22. Sensible (H) and latent (L) heat flux densities calculated for 
Death Valley for January 18, 1977. 




tABl.K n 




ATMOSlMiRHlC STABILITY PARAMRTKRS 



Condi ti on 

Parametei' 

Expression** 

1. 

Uitstable (C 0) 






0.74(1-90"^^^ 




0.74(1-9C)"^^^ 

2. 

Neutral (c " 0) 


1,0 




0.74 



♦q 

0.74 

3. 

Stable (C> 0) 


1 f 4.7C 



% 

0.74 f 4.7c 



♦q 

0.74 f 4.7c 


IkiBlng&r (1971) 

Ibly, the profile method should ftlve a closer representation of the surface 
turbulent flux densities at a point since It utlllKes not one, but many 
measurements over the thickness of the surface layer. Hiis Is particularly 
relevant at nlRht when the thickness of the surface layer is often no more 
than the height of the profile system ('^ 8 m). Furthermore, profile calcula 
tlons have been found to correspond quite well with direct measurements of 
turbiilent fluxes by the eddy correlation method (Miyake et al . , 1970), 
Practically, the bulk aerodynamic method is more convenient than the profile 
method. The latter, unlike the former, requires data derived from a complex 
array of observational equipment that Is not always easy to calibrate, 
transport, or maintain. Figure 23 shows sensible heat flux densities 
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H (PROFILE) 





calculated from both the bulk aerodynamic and the profile methods for a common 
data set collected In July, 1978 at Stonewall playa near Coldfield, Nevada. 

The bulk calculations have been made from observations at 2 m. Although the 
overall correlation between the two methods is good, there are some notable 
differences. Nighttime values (positive) of H(bulk) are larger than 
H(proflle) values. Daytime values (negative) of |H(bulk)| greater than about 
200 W/m^ show the same result. 

Tliese discrepancies occur primarily because In the bulk method the 
drag coefflent, C^, 1) Is not a function of stability, 2) Is equal for both 
momentum and heat transfer, and 3) pertains strictly to raodarate or strong 
wind conditions (cf. H' rke, 19A3; Kossby and Montgomery, 1935). The latter 
point is particularly relevant to the nighttime differences because nighttime 
winds are often light in strength, and one would expect the bulk method to 
overestimate H. This is not a serious discrepancy because nighttime flux 
densities tend to be small. On the other hand, daytime flux densities calcu- 
lated from the two methods arc relatively large and differ by as much as a 
factor of two. The net effect will be an overeat Ifiwitlon of the value of the 
thermal Inertia calculated from the model. 

Different kinds of meteorological data can be used in the model. 

For example, the National Weather Service (NWS) has a large network of weather 
stations in the continental Uulted States that routinely record weather 
observations at least once an hour. This is an attractive source of data 
since an HCMM pass is usually over at least one NWS station. Unfortunately, 
the NWS data may not be entirely representative of the weather everywhere In 
an area. Figure 2A shows two curves representing the diurnal variation of 
sensible heat flux density at two sites in the Walker Lane area of Nevada. 

The solid line depicts H calculated from NWS hourly observations in June, 1976 
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SENSIBLE HEAT aUX DENSITY (W/m'^) 


ORlQiNAL PA3E 13 
OF POOR QUALITY 


JUNE 14-15, 1976 








TONOPAH AiRPORT 

STONEWALL PLAYA 


TIME (HRS) 

Figure 24. Sensible, heat flux densities calculated from bulk 
aerodynamic method with data sets from two sites; 
Tonopch airport and Stonewall Playa, June 14-15, 1976, 




at the Tonopah airport and the dashed line represents calculations derived 
from In situ obserNatlons for the same time at Stonewall playa^ a dry lake 
about 60 km south of tonopah. Both calculations were tnade from the bulk 
aerodynamic formulations. Note several major differences. AlL'j.ough both 
curves have the same magnitude at their daytime peaks, they are not conform" 
able In general. Nighttime values of H are the roost disparate* Tlie H values 
from the NWS data are about twice as large as the M values from the In situ 
data. This pair of curves Illustrates the problem encountered when applying 
microscale measurements to mesoscale models. 

Our calculation of thermal Inertia values from the HCMM data 
Included only the effects of sensible heat transfer In the JPL model and not 
the effects of latent heat transfer because most of the areas of interest are 
seml-arld. The sensible heat flux densities have been calculated by the bulk 
aerodynamic formulation In all cases. This formulation Is sufficiently 
accurate for the calculation of thermal Inertia values from the JPL model In 
view of the Intrinsic difficulties In modeling mesoscale meteorological 
fields, and the practical difficulties of eroploylng more complicated methods 
such as the profile or eddy-correlation techniques. 

The Effect of Soil Moisture on Uterroal Inertia 

Thermal Inertia varies with the amount of water present In the 
soil. Ihls Is because the thermal Inertia of moist soil Is typically about 
0.040 cal cro“^s”^^^°C'"^ which Is comparable to thermal Inertia values for many 
rocks (cf. Kahle, 1980). As soil water content varies, the thermal properties 
of the soli also vary. The difference between the daytime maximum and the 
nighttime minimum In temperature is greater In the near surface layers of dry 
soils than It Is In wet ones. The consequences of water transport on the soil 
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thernwl risglra^ can be studied . 5 the use of nn {ippropriate mathematical 
model* 

With support from the AgRISTARS program, the JRh thermal model has 
been modified to describe heat and moisture transport In soils. Tills version 
of the model has not been used in analysis of the HCMM data for this report; 
however, it nuiy be employed in future studies. A few results friim this model 
are shown here because they illustrate the magnitude of the effect of changing 
soil moisture on the thermal history of surface nuiterlals. 

Tlie modifications to the model Include an additional term (Q) in the 
heat conduction equation to describe heat transfer by vapor diffusion under 
moisture gradients, an additional partial-differentlal equation to describe 
moisture changes, a set of initial conditions for the volumetric soil moisture 
profile, and boundary conditions to specify the moisture variation at the 
surface and at the lower boundary 1 m) of the model. Tlie result is a set 
of coupled, diffusion-type equations, 



<31) 


< 32 ) 


where T is soil temperature, 0 is volumetric moisture content, Z is depth, X 
is thermal conductivity, D 0 Is moisture dlffuslvity, D»ji is thermal dlffus- 
Ivity, K Is hydraulic conductivity, C is volumetric heat capacity, and Q is a 
term due to heat transfer by dlstlllative effects. The pair of equations is 
solved numerically by a finite difference technique. Complete details of the 
model are given in Njoku et al. (1980). 
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As « test, the moHel wn« used to slmMlnto tempernture and (nolsture 
profiles for a bare field near Ikikercf ield, California during May> 1978. 
Figures 2*> and 26 show soil temperature and volumetric moisture variations 
with depth and time. Tlie results of the modeling are shown In Figures 27 and 
28 for con^jarison* Ibe computed temperature and moisture curves agree In 
general with observations although there Is less conformity between the 
moisture curves. TIUs is not unexj>ected alnce soil moisture measurements are 
more difficult to make than temperatue measurements and ivlso, model results 
represent discrete levels In the soli while the observations are oi>lalned from 
soil layers of finite thickness. 

The variability of thermal inertia with changing moisture content 

can be Inferred from these results. Initially, the soil Is quite moist 
3 3 

(0 It 0.25 cm /cm ) to a depth of about 30 cm and the difference., AT, between 
the daytime maximum and nighttime minimum In temperature at the surface is 
about 24-29‘'*C. As the soil dries, which It does quite rapidly near the 
surface, AT increases until It reaches about 34^C at the end of the record. 
This amounts to a change In AT of about 20%. Tlie effect this change has on 
the thermal inertia can be estimated by considering the expression for 
apparent thermal Inertia (AP) in Its approxinuite form (see section V^G), 


AP S (l-A)/AT 


(33) 


where A Is soil albedo. Ignoring changes In A and taking the partial 
derivative of AP with respect to AT, one obtains 


3(AP) 

3(AT) 1 

AP " " AT 


(34) 
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MAY 1978 


Figure 25. Observed soil temperature var lot Ions at o six different depths on a 
bare r'ield near Bakersfield, California for the period May 17-19, 
1978. 



17 18 19 

MAY 1978 


Figure 26. Observed soil moisture content at six different depths on a bare 
field near Bakersfield, California for the period May 17-19, 1978. 


92 








which implies thrtt rt 20% Incrense In AT ciiuBes n 20% decrease in AP. In 
reality, A also decreaBes with Burfiice wetness, and this additional effect on 
AP iB given by 



For,' the range of 0 shown in both Figures 26 and 28, this gives an increase in 
A of about 10 to 30% as the soil dries according to Janza (1975), and Implies 
a further decrease in AP by about the fiame magnitude. Therefore, in this 
case, the thermal inertia of the soil can be reduced by nearly a half in Just 
a few days due to the effects of evaporation and internal drainage* 

Spatial and temporal variat-ons in soil moisture can be ascertained 
frora HCMM thermal Images. For example, Figures 29 to 39 show how viiese 
changes appear at different scales for the Death Valley site. The first set, 
Figures 29 to 31, are full-frame, daytime IR images of a part of the western 
United States. Death Valley is located south of the letters "DV” on the 
images of May lA, 1978, August 18, 1978, and September 19, 1978. Tlte 
elongated feature is the valley floor, the ring-shaped feature at the north 
end is Cottonball Basin designated with a "C*, south of this is Middle Basin 
shown with an "M", and further south the Badwater area is designated with er 
”B". The dark (cooler) areas in the valley coincide with the moist part of 
the saltpan. Seasonal and spatial variations in the area's ground moisture 
are evident from one frame to the next. In May, the valley is barely 
perceptible, particularly the Cottonball Basin area; in August, it is more 
distinct and Cottonball Basin is seen as a faint quasi-circular area; in 
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photograph 



KiKuro 29. Fiill-franu- IICMM day IK Image of part of tlie western United 
States for May 14, 1978, Death Valley Is located south of 
the "DV." Test site locations within the valley are Indi- 
cated by C (Cottonhall lias In), M (Middle liaslri), and li 
(lladwater) . 
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Figure 30. Full-frame HCMM day IK Image of part of the western United States for 
August 18, 1978. Death Valley test sites are indicated. 
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^>op( omiHM', du' throo lU'OiiH (('.» M iind H) liro quUt' vlwihlts parHoutfirly tho 

HiUlW.it I* iUVll. 

K1ku»*''» V*’ to an* auh.iroaa ol‘ the UOMM data that Itavo boon 
pnu'OHHod ill vIHh to dlMpUiy ttu'rmal Inortiii of Poath Valloy on Ji larp,or HoaUv. 
Thia piHHuoiaiap, haa boon doao with tho dry voraioa of ttu* dPh tiionmil im>doI. 
.Siam* of thoiu* imanoH iiro inoomploto bocamio of .a lack of ovorlappfnR HOMM data 
on amH'**aalvo panat'a; nmiotho I ohh, tlio valloy fa d laoomlblo on all fmanoa. 

Kour diiya iiro ahowas M;iy HI, l‘)7H, duly 22, 197H, Aoppiat IH, 197H, and 
April 4, iki tho May intapo, (kufmvball Hjiain la Indfcafod by tho rlap- 

ahiipod fo.it m*o noar tho oontor of tho imjipo, tho Mlddlo Haain i» tho Itpht 
aroa at tho ond of tho fllawont protrudfnp howo diatanoo from tho bottom of 
tfio rinp, and Hadwator Hiiain la looiitod by tlio hrlpiit, quaHl’-trlangular 
foiitoro an oqiml cUatiinoo bolow tho Middlo Haaiiit 'IboHo aroaa can bo 
Hlmllarly idontlfiod on tho othor fmapoa, Haoh of thoao aroaa imdorpooK 
ohanpoa in alr.o and brfphtnoHH, vialblo on tho aoparato Iwipoa, and oach aroa 
la ii ropfon of afpnffloant molatviro variability. Inoroaaod motaturo oontont 
roaul' a in hiphor thormiil fnortia, honot> wottor an*aa appoar Ughtor In tho 
imapoa. I’arta of all thoao baaina aro anb.1oot to aoaaonal flooding, and 
Hadwator Haaln oontatna rogiona whoro tho ground ia noarly alwaya aaturatod. 

U ia unfortunato that tho do lay In dlatrlbntion of IICMM data mado. it 
Impoaaiblo to idontify aroaa of algnifiofint ehango on imagoa, and thon 
undortako fiold inoasnromonta in koy aroaa during aubaoquont aatollito 
ovorpaaaoa, aa had boon originally piannod* 

Tho dotail in tho HCMM imagoa ia good onougli to idontify rogiona 1 
ohangoa in ground moiaturo but ia not good onough to dolinoato amal! foatnroa 
of apoeial Inloroat. Airorfift data aro more uaofui for Mvia taak. Figiiros 3b 
through 39 alu>w vtathlo and thormai inortia Imagoa of Poiith Valloy obtained 
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lir.MM thiMm.ll liuMt i.i im.i>',i' ot lU'.U h V.ilU'v tor Au>;ust 18. H'78. 
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WHITE photograph 


Figuro 3'j. IICMM llionnal inert ia lma^;e of Deatli Valley for April 4, 1479. 
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from aircraft data. Hu' InwRos portray the northern and southern halves of 
Death Valley. 

The visible image of the northern half (Figure 3f>) Is distinguished 
by Cottonbali Basin. Tlie Furnace Creek area Is located In the lower right- 
hand corner and the building complex and airstrip can be scon In dark, gray 
tones. Cottonbali Ikisin is a salt-encrusted portion of the valley floor and 
is subject to flooding by the rise of ground water during wet periods and by 
runoff. Tlie most notable feature in the Image is the bright area north of 
Furnace Creek on the east side of the Ikisln. Tliis feature is a very moist 
area, which also appears on the thernuil inertia inwige (Figure 38) as a bright 
spot (indicating high thermtil inertia values). On the visible inwige of the 
southern half (Figure 37), the most conspicuous feature is the dark quasl- 
trlangnlar area known as the Devil's Coif Course. It is an elevated massive 
rock salt foriwitlon which Is nearly saturated, and the surface is extremely 
rough. Tbe therimhl inertia image (Figure 39) shows this as an area of 
relatively high thermal Inertia as would be expected of a wet area. An 
vinusual feature does appear in the Middle Basin at the top, left-hand side of 
both images. Tlrls is a strip of salt-encrusted sand and silt shaped like the 
letter Y. Hie bottom leg and the right-hand fork of the. Y show a uniform 
brightness on the visible image. In the thermal Inertia image, however, only 
Che bottom half of the leg indicates high thermal inertia. Tlie difference 
between the two representations indicates the effectiveness of thermal inertia 
techniques for differentiating materials with similar surface properties but 
apparently different bulk properties. Field checking will be required to 
ascertain the nature of features in both the northern and southern halves of 
the valley. 
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m2s visible image of the souLhern 
half of Death Valley obtained from 


Figure 37 


NP-3A flight on March 11, IV77. 
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Figure 38. Tliormal Inoriia image of t lie northern lialf of 

Death Valley proLluced from data obtained during 
the NF-3A fliglit on Marcli 1 1 , 19 77. 
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Kl>*iiro 34. Tl\t'niuil liuTt i;i lm;u’r of t ho sout horii li;i 1 I of Doatli 
Viilloy proilui'od from data ohtainod daring; tho NT- )A 
f I l^ht on March 11, 1477. 


107 


M( uU *l St lu Uu is 

Tlu*rmfi 1 tnerlln c« Iculat ions dopond upon winy fnctorsf *010 model 
roHults are sonfiitlvo to varlationR in «uch pararoeters a« day-nlRlir temper'-* 
titure dlfferenee (AT), albedo, topoRraphy, sky radlanee, cloudlneaa, vegota- 
ttve cover, and ground-'level meteorology* Tliese and other related topics are 
diacuBRed in thin aectton. In particular. Table lA presents the results of a 
study to determine the magnitude of these effects on thermal inertia calcula- 
tions based on data for the Cioldfield, Nevada and Tlsgah, California test 
sites. Goldfield is about 3 ^ of latitude north of Tisgah Crater. Tlie results 
are 80iiH,‘What dependent upon the site because of the difference in Initial and 
boundary conditions in the calculations. Changes in T arc given in 
cal 

The simulations for Goldfield and Tisgah have been carried out with 
vinta from the October, 197 ^? experiments. The effects repreBent perturbations 
upon a standard value of T for each site. As noted above, the initial and 
boundary conditions — which include factors such as local meteorological 
conditions, albedo, and topography — are different for each site and tliere- 
foro the same amount of change in a source parameter can produce effects of 
different magnitudes. For example, a l^C error In measured day-night temper- 
ature difference (AT) produces an effect about four times greater at Goldfield 
than at Plsgah. Nonetheless, these arc typical results and are Indicative of 
the relative strengths of these effects in introducing errors into the 
nnalysln. 

Most of the sources of errors listed in the table are due to the way 
energy flux densities arc partitioned at the upper boundary, or surface. For 
Instance, the albedo, ground slope and azimuth, atmospheric radiance (a.k.a* 
sky temperature component), and cloud cover all affect the composition of the 
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tARtR U 


PERTURBATIONS UPON THERMAL 

INERTIA MODELS 


Source 

EEFECT 



GOLDEIELD 

PISGAH 

l^C Error in Measured AT 

.012 

*003 

2®C Error In Measured AT 

.023 

.006 

8Z Error in Measured Albedo 

.oil 

.007 

A® Error in Measured Ground Slope 

.005 

.002 

A5® Error in Measured Ground Slope Azlmutlt 

.013 

.005 

20°C Error in Estimated Sky Temperature 

.001 

.003 

Effect of 20% Vegetation Cover 

.014 

.006 

Effect of Morning Cloudiness 

.001 

.002 

Effect of Noon Cloudiness 

.008 

.007 

Using Ely, Nevada Weather Station Data 



Wind Speed Only 

.009 


Air Temperature Only 

.047 


Both Wind and Temperature 

.026 



net radiative flux density at the surface. Errors due to subsurface factors 
such as water transport in soil? are not considered here. Cloud cover is 
simulated by eliralnatlng the direct solar heating for a short time during 
early (morning) and late (afternoon) phases in the diurnal cycle. Only the 
diffuse component of radiation incident at the surface is included in the 
calculations. The cumulative effect of errors in the radiative factors is 
considerable both at Goldfield (0.053 cal cm’*^s''^^2 °C"^) and at Plsgah 
(0.026 cal cm"^s“^'^^ °C"b» and they rival the sum of the remaining effects. 
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Veg^trttlvc. cover 1 b Important because It not only nffecta the areal albedo ], 
but also the near surface micro -cllimite and tarbulent heat transfer. However, 
In our analysis, we have only conBlderod the effects of temperature modiflca*" 
tlon by assuming that the ground temperature of the vegetated areas was e<iual 
to the air temperature. 'The final entry under the Coldfield column Indicates 
the Important role that sensible heat flux ploys In redistributing energy at 
the surface. Sensible heat flux Is proportional to the product of the wind 
and the air-ground temperature difference. We conclude that although errors 
In both wind and temperature seriously affect the calculation of P, an 
erroneous measurement In the air temperature is particularly damaging. Over 
natural surfaces this can be a major problem because air temperatureB often 
vary substantially on a regional basis, as discussed earlier. Errors In the 
measured temperature also create large uncertainties in calculating P. Tb’ise 
errors are multipUcatlve} l.e. a 2^C error in AT produces about twice as much 
error as a I'^C error. Their size is brought into perspective if you consider 
the inherent error of 0.5^C In an IR radiometer like the PRT-5. 

The results In this table Indicate that errors can be Introduced 
into the thermal Inertia modeling through measurement (e.g. temperature, 
albedo, topography) or parameterization (e.g. sensible heat flux, cloud cover) 
and that these enors can account for spurious contributions to P of signifi- 
cant magnitude. Despite the effect of these errors, the model gives results 
that are accurate enough to discriminate among many different kinds of 
geologic materials. Improvements In the model's accuracy will depend upon 
further research In technical and theoretical areas related to measurement and 
paramecerlzatlon techniques. 
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Numerical Vaiuea of ThennaT Inertia 

Table 15 la a list cl thermal Inertia values and their standard 
deviations for selected small areas of fairly uniform composition at the four 
test sites. All the data pertain to model calculatons of P obtained from the 
July 6, 1978 data sets. Also listed are thermal Inertia values for some of 
these materials as published In the Manual of Remote Sensing (Reeves, 1975) In 
the column headed MRS. In an absolute sense, our values of P do not always 
agree closely with the limited published values for the materials. In general 
our results are larger In magnitude. However, in a relative sense within a 
given site at a given time and v;ithln the limits of error, the values of P In 
Table 15 are useful In discriminating among many rock and soil types. For 
example, at Coldfield, altered and unaltered materials are readily distin- 
guishable from each other, but almost all values appear to be too high. At 
Plsgah, different kinds of basalt can be separated one from another, while at 
Death ValJey, dolomite Is separable from other materials. It is Interesting 
to note the wide range of values of thermal Inertia among the sandstone units 
at San Rafael. These differences appear to be real and are ascribed primarily 
to the different weathering properties and resultant soli development on the 
various units. 

The fact that our thermal Inertia values are In general larger than 
published values can be ascribed to our not yet Included atmospheric correc- 
tions. Because of the surface temperature dependence of these corrections, 
the warmer daytime temperatures will be more affected than the cooler night- 
time temperatures. It is apparent from Table 9 that the measured AT, with no 
atmospheric corrections, will be smaller than the actual AT. For Instance, an 
actual daytime surface temperature of 320°K and nighttime temperature of 290°K 
(Al “ 30°) would result In a satellite measured AT of only 24°. This smaller 
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TABLE 15 


THERMAL INERTIA VALUES AND THEIR STANDARD DEVIATIONS AT FOUR TEST SITES 


P 0 


MRS Value 


7/6 Goldfield 


7/6 Plsgah 


7/6 Death 
Valley 


San Rafael 


1 Playa 

30.95 

5.00 

A2 (Clay soll-molst) 

2 Basalt 

37. 8A 

5.09 

53 

3 Slllclfied 

5A.67 

2.58 


Volcanlcs 




^ Argllllcally 

35.50 

1.06 


Altered 




Volcanlcs 




5 Unaltered 

AA.31 

1.86 


Volcanlcs 




1 Playa 

53.96 

3.A8 

A2 (Clay soll-molst) 

2 Aa 

81. A7 

2.63 


3 Basalt 

76.73 

1.77 

53 

4 Basalt 

7A.62 

3.65 

53 

5 Blotlte Quartz 

77.33 

3. A3 


Monzonlte 




6 Alluvium 

61.22 

1.13 


7 Daclte 

55.87 

1.78 


8 Granite 

7A. 13 

1.02 

52 

1 Dolomite, Qtz, LS 

10A.93 

6.17 

75, 7A, A5 

2 Dolomite, Qtz, LS 

107. OA 

6.15 

75, 7A, A5 

3 Quartz Monzonlte 

99. 7A 

5.8A 


A Alluvlun, Pleistocene 

78.87 

2.18 


5 Q Lake 

60.62 

2. AO 


6 Alluvium, Recent 

86.01 

1.79 


7 Sand 

61.26 

1.72 

2A (sandy soil) 

8 Basalt 

82.05 

A. 07 


1 Navajo SS 

72.76 

A. 20 

5A 

2 Entrada SS (soil) 

9.22 

2.30 

5A 

3 Kalbab LS 

70. 6A 

8.86 

A5 

A Moenkopl Silt 

A9.56 

A. 22 


5 Cedar Mtn Fm 

31.56 

1.95 


6 Carmel Fm 

18.39 

1.89 


7 Mancos Shale 

20. A9 

3.A6 

3A 

8 Perron SS 

27.33 

2.27 

5A 

9 Salt Wash SS 

A0.2A 

5.37 

5A 

10 IXines 

27.02 

2.06 

2A (sandy soil) 
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AT would result In erroneously larger values of calculated thermal Inertia* 

As seen In Table 14, even a 1° error In AT can result In an error In thermal 
inertia of from .003 to .01 cal cm”^s“^'^^°C” ^ or more. This clearly demon- 
strates the Importance of the atmospheric corrections which we will be In- 
corporating Into these calculations. 

Large variations in calculated values of thermal Inertia occur at 
all sites, as a function of time. Table 16 shows the average values of 
thermal Inertia (P) and the standard deviations ( o) for the HCMM overflights 
of 1978. The top set of numbers Is calculated from the JPL tLermal model and 
the bottom set from apparent thermal Inertia calculations (cf. section G for a 
discussion of apparent thermal Inetrtla). Values of P refer to the entire 
Image and are obtained from calculations of the frequency distribution of P 
over the scene. The bottom three rows In the table give the values for all 
days of the range of P, the average value of the ratio of 0 to P, and the 
average value of P, respectively. Goldfield displays the greatest amount of 
variability for both cases. This list gives some Idea of seasonal changes 
in P. Both Death Valley and Pisgah data show increases or decreases in P with 
time, and there is also some indication of an Increasing trend in P in the 
Goldfield data although this Is not explicitly clear because of the variabil- 
ity in the results. The model results show the changes in P at Goldfield from 
one date to the next are almost as great as the range in P over the entire 
sample. 

Note that the model values of P are larger than the apparent therme.1 
Inertia values in all cases although the general form of the variations in 
P is similar throughout. 

Seasonal changes in the calculated value of thermal inertia can be 
attributed to a number of factors, e.g., variations in meteorological condi- 
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TABLE 16 


AVERAGE VALUES (P) AND STANDARD DEVIATIONS (a) 
OF THERMAL INERTIA (x 1000) FOR 1978 HCMM FLIGHTS 


MODEL CALCULATIONS 


Date 

Goldfield 

Death Valley 

Plsgah 

May 14 

NA*** 

89.5 k 21.5 

64.7 k 10. 

May 30 

NA 

81.6 k 14.4 

65.6 k 9.1 

July 6 

50.0 k 17.6 

76.6 i: 13.4 

63.8 k 8.0 

July 17 

45.0 k 15.5 

76.7 k 13.5 

71.0 k 9.4 

July 22 

62.5 1: 15.2 

77.4 k 16.3 

69.5 k 20. 

Augu St 18 

36.5 k 14.5 

73.1 k 13.9 

72.6 k 9.1 

September 19 

67.5 i: 13,2 

NA 

84.3 i: 7.4 

Range 

17.5 

16.4 

20.5 

<o/P>* 

0.31 

0.19 

0.15 

<P>** 

52.3 

79.1 

70.2 



APPARENT 

THERMAL INERTIA 


Date 

Goldfield 

Death Valley 

Plsgah 

May 14 

NA 

53.1 ± 12.8 

38.8 J: 4.6 

May 30 

NA 

45.1 k 8.4 

38.7 k 4.2 

July 6 

39.8 k 8.4 

46.2 t 9.3 

38.2 J: 3.6 

July 17 

36.5 k 5.5 

46.8 k 8.5 

45.0 i- 5.0 

July 22 

38.9 k 9.2 

47.8 k 10.7 

44.3 k 8.1 

August 18 

35.1 i: 5.1 

47.5 Jr 8.8 

42.6 k 4.9 

September 19 

1 44.9 k 10.3 

NA 

47.4 f 5.6 

Range 

9.8 

8.0 

9.2 

<a/P> 

0.19 

0.20 

0.12 

<P> 

39.0 

47.7 

42.1 

* <o/P> 

= average o^ o/P's 



** <P> = 

average of P's 



*** NA = 

not analyzed because of 

Incomplete data 
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tlons, in ground moisture, in surface albedo, in vegetation cover, and in the 
emission and transmission properties of the atmosphere along the IR signal 
path length. Some of these factors cause true changes in the thermal inertia, 
and some only Introduce errors into the calculations. 

We have investigated the effect of two factors, local wind speed and 
atmospheric water vapor content. We first considered variations in the local 
winds. We have calculated the correlation coefficient (r) for the average 
daily wind speed and the values of P for each day at each site listed in Table 
16. For both model and apparent values of thermal inertia, no significant 
correlation (r ~0.l) exists between wind strength and P. This result is 
partly due to the fact that the wind speeds used in the analysis were measured 
at one location for each site whereas the mesoscale variability of winds over 
an entire site should be considerable. 

Second, we have compared values of P with the amount of water vapor 
present in the atmosphere. Tills analysis has only been done for the Goldfield 
data because Goldfield is the only test site located near a rawinsonde 
station. Rawinsonde data were obtained from the Mercury, Nevada National 
Weather Service upper air station which is about 160 km southeast of Gold- 
field. Tlie upper air data consist of pressure, temperatvire, and dew point 
measurements for a number of heights from near the surface to about 30 km. 

The amount of preclpltable water (W) in a column of air of unit cross-section 
can be calculated from the upper air data (Saucier, 1955) and, therefore, W is 
a measure of the amount of water vapor present in the atmosphere. W is an 
important parameter to consider because water vapor absorption is the princi- 
pal mechanism for attentuatlng the IR signal strength along its path of propa- 
gation. Nonetheless, we did not find any significant correlation (r “ 0.10) 
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between P rtnd W from the Ooldfleld and Mercury data for either the model or 
apparent thermal Inertia calculations. 

Ground moisture Is another factor of Importance when considering 

variations in P between Images of the same site. K few measurements of soli 

moisture content are available from the sites, but the areal and temporal 

coverage during the aircraft and HCMM overflights was not sufficient for 

detailed study. Furthermore, most of the measurements pertain to very dry 

3 3 

condltons - 0 < 0»05 cm /cm' - and so a statistical analysis of the data would 
not be meaningful. 

The simple correlations In the preceding text may not be appropriate 
because the many variations In P interact in an Intricate way. Soil moisture 
effects undoubtedly play a role in the variability of P, but they may be 
modified by other condltons. In complex terrain, changes In P due to soil 
moisture compete with changes engendered by warm or cold air advectlon. The 
latter effects are different than the turbulent transport of heat between the 
ground and air which Is accounted for in the model calculations of P. Those 
calculations are made under the assumption that advectlon Is unimportant. In 
reality, mesoscale wind systems often advect air of different thermal proper- 
ties from one place to another. Of particular Importance In our analysis are 
drainage winds which usually occur during the night, and are also known as 
slope or valley winds. Ihese winds frequently result In cold-air ponding on 
valley floors and In Isolated hollows In open country. The calculated thermal 
Inertia of materials that are located In depressions and are subject to cold- 
air ponding may be underestimated because the day-night temperature difference 
Is Increased above the value it would obtain In the absence of these condi- 
tions. 
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As an example, consider a circular basin of uniform material subject 
to cold-air ponding. Tlie calculated value of P at the center will be lower 
than that at the edges so that after processing, the corresponding Image will 
show an annular band of high P surrounding a central core of low P, This 
result Is complicated by the fact that water collects In depressions, and this 
has the opposite effect on P. P at the bottom, or center. Is then higher with 
respect to Its value at the drier edges, and so a thermal Inertia Image would 
show a light center enclosed by a dark ring. This Is Illustrated in Figures 
40 and 4l which are thermal Inertia Images of the Goldfield area. The upper 
right-hand corner of each Image shows a quasl-clrcular playa known ts Mud 
Lake. Ihls appears as a bright disk surrounded by a dark annular ring on the 
April, 1979 Image. This Is indicative of the latter condition discussed In 
which the outer edges of the basin are dry and the center is moist. The same 
feature Is visible on the July, 1978 image except a small dark area is super- 
imposed upon the center of the bright disk. This is circumstantial evidence 
of cold-air ponding in a moist depression. This view is supported by a review 
of wind measurements that were made on nearby Stonewall playa. Steady, 
moderate winds (more than 2 m/s) occurred In the early morning hours of April 
4, 1979 in the area and Intermittent, weak winds (less than 1 m/s) occurred 
for the same time on July 6, 1978. The latter winds usually Indicate drainage 
wind In this kind of terrain. 

G. Comparison of Apparent and Modeled Thermal Inertia Images 

Apparent thermal Inertia (AP) Is defined by Price's (1977) formula, 

AP « 2 SVC . (1-A) (36) 

Ml ^ f 
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Figure 40. Tliermal inertia image obtained from HCMM data for April 4, 1979 in 
tlie Walker I.ane area near Coldfield, Nevada. 


118 




ORIGINAL PAGC 

BLACK AND WHITE PHOTOGRAPH 


Figure 41. Tliermal inertia image obtained from HC^D•t data for .luiy b, 
1978 in the Walker Lane area near Goldfield, Nevada. 


119 


where S is the solar conetant (- 1.98 cal cm”^ V la atmospheric trans- 

mittance at visible wavelengths (a value of 0.75 is used In the calculation), 

A Is albedo, AT Is the day-night temperntu/e difference (taken between 1330 LT 
and 0230 LT), o) Is the angular speed of the earth's rotation (- 4.363 x 10 
radlans/mln) , a Is a parameter which Is proportional to the ratio of the heat 
flux density transferred by the surface to the air to that transferred Into 
the ground, and C Is a parameter which varies with solar declination (A) and 
latitude ((^) according to 

C * i [slnA sin^ arccos (- tanA tan^) 

f cosA co8^(l-tan^A tan^^)^'^^] . 


'’’able 17 shows the result of a linear least squares fit of P values calculated 
from the JPL thermal model versus AP values calculated from Equation 36. 
Meteorological data from four midlatitude NWS locations - Ely, Nevada and 
three Arizona stations, Yuma, Winslow, Tucson - have been fit to the linear 
equation. 


P(model) - mAP f b. (38) 

The last column In the table gives standard deviations ( o) for each calcu- 
lation. The seasonal values refer to solar declinations of A » 23.5° (June 
21), A = 0° (March 21, September 21), A ” -23.5° (December 21). 'fhe model 
calculations pertain to dry conditions (no soil moisture or latent heat flux 
components) for 96 cases: four stations x four seasons x three P values 

(0.02, .05, .08 cal cm~^ °C”^) x two albedo values (0.1, 0.4). Values 
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TABLE 17 


LINEAR LEAST-SQUARES FIT OF P( MODEL) TO AP 



m 

b 

0 

ALL CASES 

2.09 

-0.026 

0.012 

Ely, Nevada 

1.61 

-0.008 

0.017 

March 

3.06 

-0,055 

0.005 

June 

2.92 

-0.080 

0.008 

Sept 

3.12 

-0.048 

0.004 

Dec 

4.29 

-0.057 

0.009 

Yuma, Arlzo't.a 

2.61 

-0.044 

0.006 

March 

2.78 

-0.052 

0.004 

June 

2.69 

-0.052 

0.003 

Sept 

2.54 

-0.042 

0.003 

Dec 

2.77 

-0.043 

0.005 

Winslow, Arizona 

2.31 

-0.032 

0.012 

March 

3.02 

-0,048 

0.004 

June 

2.87 

-0.065 

0.004 

Sept 

2.98 

-0.062 

0.007 

Dec 

3.09 

-0.041 

0.005 

Tucson, Arizona 

2.64 

-0.050 

0.006 

March 

2.81 

-0.053 

0.005 

June 

2.49 

-0.048 

0.003 

Sept 

2.76 

-0.059 

0.006 

Dec 

2.91 

-0.054 

0.005 
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of AT arc computed for each case and^ along with other data, arc used to 
compute AP values which arc used with the given P values as data in the 
regression analysis* Different values of m, b, and 0 occur for similar 
seasons (e.g* March, September where 6-0) because the model calculations 
employed climatological data for each station, and these data differ from One 
month to the next. The values of m and b for all cases are 2.09 and **0.026 
respectively and the table shows that P(model) is larger than AP in all cases. 

It is Interesting to apply this result to the P(model) and AP 
calculations of the HCMM flights listed in Table 1. Tlte AP values of the HCMM 

data are also lower than the corresponding P(model) values. If the HCMM data 

for AP are used in the formula 

P(modcl) « 2.09 AP ~ 0.026 (39) 

then the values of P(model) calculated from this formula agree, for the most 
part, within 5 to 15% of the P(model) values listed for the HCMM data. 

Two Versions of the modeled Images of P are used*, one with and one 

without corrections for topography. The topographic corrections are obtain ,d 
from NCIC data which include Information on elevation in digital format. 
Figures A2, A3, and AA show topographic images derived from the NCIC data of 
the Goldfield, Nevada area. These are respectively a shaded relief map, a 
slope azimuth map, and a slope orientation map. For reference, the Goldfield 
mining district is in the upper right of each image. Stonewall Mountain is in 
the lower right corner of each and Clayton Valley is in the left central 
portion. 

The next set of figures, A5 through A7, show HCMM data for July 6, 
1978. In order, they are apparent P, modeled P without topographic correc- 
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Figure 42. 


Shaded relief map 


of the Walker Lane, 


Nevada area. 
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Flyure 43. Slope azimuth map of the Walker Lane, Nevada aiea. 
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P’igure 44. Slope orientation .nap of the Walker l.ane, Nevada area. 
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tlons, and modeled P with topographic corrections. The large, white area In 
Clayton Valley Is a body of standing water. There are no significant differ- 
ences between the modeled Images. Evidently, the NCIC scale (~ 60 m) for the 
topography Is too fine-grained to be useful on the scale (~ 500 m) of the HCMM 
data. Furthermore, the apparent P Image also displays little variation from 
the other two Images. This also holds true for the other HCMM data products 
tikat were processed but not shown here. 'Rie one exception noted Is at Mud 
Lake which Is seen as a dark or gray disk-like feature In the upper right on 
each image. The dark area around the perimeter of the disk Is more extensive 
on the apparent P Image than on the modeled one. The reason for this differ- 
ence Is not clear but it may be related to the ground moisture distribution. 

This study Indicates the benefits of constructing apparent thermal 
inertia images. The quality of the Images Is comparable to that attained by 
the more sophisticated model products without the extra labor and cost 
involved. This advantage is restricted to image Interpretation in a relative 
sense, l.e. discriminating among geologic materials with different thermal 
inertia values. In an absolute sense, neither the apparent nor model thermal 
Inertia images faithfully reproduce exact thermal inertia values of the 
geologic units depicted. However, the model products are more likely to 
approach the actual values since model calculations make provision for some of 
the factors (e.g. atmospheric effects, ground-level meteorological conditions) 
that cause discrepancies between calculated and actual values. This does not 
minimize the Importance of the apparent thermal inertia products. For many 
geologic applications, the vise of thermal Inertia images that are valid in a 
relative sense is quite appropriate, and in this context, apparent thermal 
Inertia products are considerably useful. 
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l’i^>urc' 43. Appari’nt Ihonnal iiu>rt ia image ohLained 1 rom HC'.MM data for Walker I.ane 
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Figure 46. Calculated thermal inertia image obtained from HCMM data for 
Walker Lane on duly 6, 1978. Topograpliic corrections not 
inc luded . 
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F4>uru A7. Calculated thermal LnerLia image obtained t'rom IICMM data for 

Walker Lane on July h, 1978. Topographic corrections Included. 
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VI. GEOLOGICAL INTERPRETATION 


A» Walker Lane 

1. Geologic Description 

The Walker lane test site Is In Esmeralda County, southwestern 
Nevada, (Figure 48), In the Basin and Range province* The overall topography 
Is dominated by arcuate ranges and Intervening valleys. The general relief 
between valley bottoms and the crests of adjacent ranges Is on the order of 
1000 m. 

The climate Is arid to seml-arld with an average annual precip- 
itation of about 10 cm. Vegetation on the lower slopes and valleys consists 
of sagebrush and other desert plants; average cover Is less than 15%, At 

tm 

higher elevations plnon or juniper predominate. 

Rocks exposed In the area Include sedimentary, volcanic, and 
Intrusive rocks. Sedimentary rocks consist of dolomite, limestone, shale, 
valley alluvial fill derived from the adjacent mountains, and playa deposits. 
Plutonic rocks are of Mesozoic and Tertiary age, and are dominantly quartz 
monzonltlc In composition. Volcanic rocks Include welded and non-welded ash 
flow tuffs, lava flows and volcanic breccia, ranging In composition from 
rhyolitic to basaltic, with quartz latlte composition being the most common 
(Albers and Stewart, 1972) (Figure 49, geologic map)* 

The Walker Lane test site area Is famous for Its precious metal 
deposits. The Goldfield and Tonopah gold and silver deposits have produced 
major tonnages of ore since 1890. Hydrothermal alteration associated with the 
Goldfield deposit Includes nillclf Ication and argillization. To the south, 
the Cuprite mining district is characterized by slliclf ication and opaliza- 
tion. 
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Figure 48. Location map of Walker Lane, Plsgah and Death Valley sites 
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2* HG^tM Images 

An HCMM day visible Image for the Walker lane test site is 
shown as Figure S0> (Tills and all subsequent HCMM scenes lor Walker Lane, 
Death Valley, and Plsgah are from the 6-7 July 1978 data set.) The very 
bright areas are playas and a few clouds (in the northwest corner). The 
darkest areas are basalts and other volcanic rocks. At high elevations, 
vegetation cover also appears dark. 

The day IR Image (Figure 51) has Clayton Valley playa and the 
small clouds displayed as the coldest areas In the scene. Tliere Is probably 
standing water on the playa, as the other playas In the scene appear much 
warmer. Other cold areas Include Stonewall Mountain, and other high eleva- 
tions. Both the altitude and vegetation cover causr* the cooler temperatures. 
Alluvial valleys appear very warm, with low to moderate elevation bedrock 
regions appearing cooler. Very little. If any, lithologic Information can be 
interpreted from this presentation. 

The tonal patterns on the night IR Image (Figure 52) correspond 
closely to the pattern of outcrop areas vs. alluvial areas, with the former 
generally colder than the latter. Clayton Valley has a very w^rm area In the 
north (water) and a cold area to the south (dry). 

The thermal Inertia image (Figure 53) has more geologic Inform- 
ation than the previous three data types. The highest thermal Inertia areas 
(white) correspond to outcrops of sedimentary rocks of dolomltlc, quartzltlc, 
and limestone composition. Low thermal Inertia areas Include alluvial fans 
and dry playa deposits. The area around Goldfield Is mixed; a high thermal 
Inertia area corresponds to an area of basalt and slllclflcatlon; moderate 
thermal inertia areas are rhyolites and andesites; low thermal Inertia areas 
ace alluvium. 
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Figure 50. HCMM day visible image of Walker Lane site. Image size Is 
65 X 85 km. 
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Figure 52 


.ICMM niglit IK linage of Walkor Lane slLe. linage size 
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Ki^ure 53. MCMM tliernuiL ima^;f of Walker Lane site. Image size Is b5 x 85 km. 
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A far more affective display of the HCMM data Is shown In 
Figure 5A. Tlie color composite was produced by combining the visible, day IR, 
and night IR Images, displayed as green, red, and blue, respectively. The 
visible and night XR Images were also complemented. Complementing the data Is 
the same as changing a photographic positive to a negative. BlacR areas 
become white, white areas become black. The algorithm used Is: 

DN - 255 " DN, 
o 1 

where DN^ Is the digital number or brightness value of the new picture, Is 
the value of the original picture. The brightness values range from Q (black) 
to 255 (white). Areas which have low albedo, are cold at night, or are warm 
in the day, appear brightly colored In the composite picture. An Interpret” 
atlon map produced from this picture Is shown in Figure 55. Table 18 Is a 
description of the symbols used for the Interpretation. 

The playas are dark blue (^1) as are the altered rocks at 
Goldfield (fl^l6) and the altered rocks at Cuprite to the south. At Goldfield, 
alteration consists of both slllclf Icatlon along faults, and more widespread 
argllllzatlon. It Is the latter that probably stands out. The similar appear- 
an'te of playa material Is due to their similar albedos and thermal properties 
(low thermal Inertia). 

Basalts are yellow or green; other volcanic rocks are yellow- 
green. However, some outcrops of alluvium and phyllltlc siltstone (//14) also 
appear the same color. Red to pink areas are alluvial valleys. There are a 
great many areas where the different lithologies cannot be separated (H5), 
These appear In various colors. Overall, some lithologic separations could be 
made; more Information Is displayed than on any of the separate black-and- 
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Figure 5^. IICMM color composite of Walker Lane. Day IR, night IR, and 

visible images displayed as red, blue, and green respectively. 
Visible and night IR images were complemented. 
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INTERPRETATION OF HCMM COLOR COMPOSITE: DAY VISIBLE + DAY IR + 
NIGHT IR. FOR SYMBOLS, SEE ACCOMPANYING CHART 


Figure 55. Interpretation map produced from Figure 54. 
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TABLE 18 


DESCRIPTION OF SYMBOLS FOR WALKER LANE HCMM INTERPRETATION 


Unit 

GcoIorIc 
Map Symbol, 

Material 

1 

QP 

Play a 

2 

Qa , QToa 

Alluvium 

3 

QTs 

Ijake deposits 

4 

QTb, Tba 

Basalt 

5 

Qta, Ta3, Ta2 

Andesite 

6 

Tr3, Tr2 

Rhyolite 

7 

Tt3, Tt2 

Silicic ash flow tuff 

8 

Ts3 

Tuff sediments 

9 

Tri 

Rhyolite Intnislve 

10 

Tml 

Mafic intrusive 

11 

T(!?r, Jgr, KJd, Mzgr 

Granite 

12 

OC 

Limestone, dolomite, shale 

13 

OS, €t 

Sliale, limestone, chert 

14 

6Zs , Zw 

Phyllitlc siltstone 

15 


Mixed 

16 


Alteration at Goldfield 
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white timiges. However, other work using lands.'it data was far superior for 
lithologic nwpplng In this area (Rowan et al» , 1974). 

3. Aircraft Images 

High resolution aircraft data were acquired over the Goldfield 
and Cuprite mining areas. A thermal Inertia image for the Goldfield site 
appears as Figure 56. Hie same image with the mapped alteration zones super- 
imposed (Ashley and Keith, 1976) appears as Figure 57. The alteration ssoivcs 
Include areas of high thermal Inertia (slllclfied) and low thermal inertia 
(argllltzed). Other high thermal inertia areas are unaltered volcanic rocks 
(basalts, andesites); other low thernuil Inertia areas are unaltered rocks and 
alluvium. Unique Identification of the alteration zone is not possible from 
the thermal inertia data. Hie circular distribution of high thermal inertia 
material, however, would alert the geologist to the possible 'Existence of a 
structural feature. In this area, this feature reflects control of the alter- 
ation by a caldera boundary. Alteration fluids were channeled along fractures 
at the boundary, producing the circular pattern. 

Tlie geology and alteration maps for the Cuprite dlstict to the 
south of Goldfield are shown in Figures 58 and 59. An aircraft thermal 
inertia image of the area is shown in Figure 60. The slllclfied area in the 
eastern half of the district appears as a region of high thermal Inertia on 
the image. On the western half of the image, high thermal Inertia areas 
correspond to opalized and sillclfled rocks. The playa on the right edge of 
the image appears dark (low thermal inertia). T!ie appearance of the slllcl- 
fied and opalized rocks is consistent with the appearance of similar rocks in 
the Goldfield aircraft thermal inertia image. 

In order to better evaluate the contribution of thermal data 
for rock type/alteration discrimination, the thermal inertia data were regls- 


142 



BLACK AND WHITE PHOTOGRAPH 



Kij;urt' '3(i. rivcrm.il liu'rl i.i ima>’,c' of Colcifiold, Novada from aircrafl data. 
lma>;o slzo Is 15 x 15 l<Ji\ . 


1A3 



B'^XK 


A* 


PAGE 

:d white photograph 





Mining; District, Nevada 



. J ' 





,y.'. 









ORIGINAL PAGE 
COLOR PHOTOGRAPH 


taieiMU Mca 




PAGC 

PHOrOuRAPH 





Figure 60. Thermal inertia image of Cuprite Mining District, Nevada 
from aircraft data. Image size is 8 x 12 km. 
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tered to an earlier aircraft data act, obtained with NASA's 24-channel scan- 
ner. This data set Included data In the 0.5 to 2.5 |in region. A color ratio 
composite picture produced from this visible and near-infrared data (Figure 
61) allowed discrimination of the three mapped alteration types (Abrams 
al. , 1977). The visible and near-infrared data were combined with the thermal 
Inertia data, and processed using a principal components analysis (Blodget et 
al. , 1978). A color conf>o8lte produced from the analysis Is shown In Figure 
62. Slllclfled rocks appear yellow, argllllzed roci:s deep blue, and opallzed 
rocks green. No additional alteration Information was produced. However, 
after Inclusion of the thermal data, at le.lst three separate units can be 
discriminated In the unaltered volcanic rocks In the north-central part of the 
scene (brown on the geologic map. Figure 58), which are not shown on the 
geologic map or on the color composite produced from the visible and near- 
infrared data above. Field work Is necessary to Identify the causes of the 
image distinctions. It la probable that the image processing brings out sub- 
units of the volcanic rocks which either were not recognized or were combined 
on the map as a single unit. 

B. Death Valley 

1. Geologic Description 

Death Valley Is In southeastern California (Figure 48) at the 
south edge of the Great Basin. The valley trends north-south between block- 
faulted mountains. In the main part of the valley, the floor Is a flat playa, 
crusted with salts, and covering more than 500 square kilometers. The moun- 
tains bordering the valley rise to more than 3400 m only 7 km from the edge of 
the salt pan, 80 m below sea level, > 

Death Valley Is a desert area, receiving an average of 4.2 cm 
of rainfall annually. Below about 1500 m, Death Valley is part of the lower 
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Cuprite Mining District, Ncvadft . Color ratio composite from 
visible and near- inf rared aircraft mu 1 1 ispect ra 1 scanner data. 
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Figure 62. Cuprite Mining District, Nevada. Principal components composite 
visible, near-infrared, and thermal inertia aircraft data. 



Sonoran j5onc» characfcerired by creOBOte biish* Vegetation 1* ecanty or abeent, 
except along the valley floor where aprlng-fed areaa locally support nesquite^ 
salt bush» and plckleweed. Overall, vegetation density la no more than 5%. 
Above 1500 m plnon- juniper become prevalent • 

Rocks exposed In the Death Valley area range In age from Pre- 
Cambrian to Quaternary* They Include metasediments; dolomites, limestones, 
sandstones, shales, quartzites of marine origin; granitic Intruslves; volcanic 
rocks of basaltic, andesitic, and pyroclastic composition; continental sedl** 
raents and fanglome rates; evaporlte deposits of varying saline and gypsiferous 
composltloa; and sand dunes. 

Structural features are dominated by Quaternary normal Basin 
and Range faults which formed the trough that is Death Valley, and northwest 
trending right-lateral/normal faults. Older thrust faults and normal faults 
in the mountain ranges produced a structurally complex exposure pattern of the 
PrN'icambrlan to Mesozoic rocks (Hunt and Mabey, 1966) (Geologic map. Figure 
63). 

2. HCMM Images 

An HCMM thermal Inertia Image for the Death Valley test site Is 
shown In Figure 64* Bright areas have high thermal Inertia and dark areas 
have low thermal Inertia. The brightest features In the scene correspond to 
areas underlain by dolomite, limestone, quartzite, or granite. These rock 
types stand out distinctly on the thermal Inertia image and are quite easy to 
map. The darkest features correspond to areas of young alluvium in the moun- 
tain valleys while the floor of Death Valley is bright to medium fray; the 
former areas are those that are probably more moist than the other areas (cf. 
Section V-D). 

The thermal inertia image also allows accurate delineation of 
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Figure 64. HCMM thermal inertia image of Death Valley area. Image size is 
110 X 126 km. 
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the bedrock -alluvium contact. This differentiation Is often difficult to make 
on visible or reflected IR Images, particularly when moderate to low spatial 
resolution prevents the use of textural Information. 

Examination and comparison of eight thermal Inertia Images 
produced from different dates of HCMM data indicates that there are no differ- 
ences In the appearance of the mountains and alluvial areas. Some change was 
observed in the playas, as previously described. 

Color composites created from the HCMM data produced the best 
overall Images for displaying lithologic differences. Figure 65 is a combin- 
ation of night IR, day IR, and visible displayed as red, blue, and green, 
respectively. The three components have been complemented, so areas that are 
cold or have low albedos are brightly colored. 

The Death Valley salt pan Is deep blue. Indicating It has high 
albedo, Is cold In the day and warm at night. This Is consistent with a 
composition of clay, silt and salt that Is also moist. A detailed interpret- 
ation of boundaries that are dellneable on the image, and the mapped geologic 
units they correspond to, Is shown In Figure 66. A description of the symbols 
Is listed In Table 19. 

Light blue areas (#13) are outcrops of limestone, dolomite and 
quartzite. These materials have medium albedos, are cold In the day and warm 
at night; these are also the areas which had the highest thermal Inertia 
values on the previous Image (Figure 65). Green areas are underlain by older 
alluvial fans. Yellow-green areas (#5, 7, 8) are underlain by fanglome rates, 
loosely consolidated volcanic rocks and basalts. Orange and red areas corres- 
pond to young alluvial deposits (cold at night, warm In the day, moderate to 
low albedo) and are the areas with the lowest thermal inertia values In Figure 
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65.. HCMM color composite of Doath Valley. Day IR, night IR, and 
visible Image.s displayi'd as blue, red, and green, respective- 
ly. The three images were also complemented. 
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Figure 66. Interpretation map of Figure 65. 
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TABLR 19 


DESCRIPTION 

OP SYMBOLS FOR 

DEATO VALLEY HCMM INTERPRETATION 

Unit 

CeoloRlc 
Map Symbol 

Material 

1 

Qs 

Sand 

2 

Qal 

Alluvium 

3 

Qst 

Sa It 

4 


Lake beds 

5 

Qc, QP, Pc,Tc 

Conglomerate, Fanglomerate 

6 

QPvb, Pvb 

Basalt 

7 

Pvr, Tvr 

Rliyollta 

8 

Pvp, Tvp 

Pyroclastic. 

9 

Tir 

Rhyolite Intrusive 

10 

Tgr, gr 

Granite 

a 

bi 

Basic Intrusive 

12 

m 

Metamorphlc 

13 

i,D,s,o,e 

Limestone, Dolomite, Quartzite 

14 

Im 

Quartzite, Sandstone, Limestone 

15 

€? 

Limestone, Shale, Quartzite 

16 

«p6 

Conglomerate, Quartzite, Granite 



Gneiss 

17 

ep6 

Granite, Gneiss 

18 

CP, CM 

Shale, Limestone 

19 


Mixed 
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The distinction between alluvial areas and bedrock is particu- 
larly well defined. In addition, some separations can be made between differ- 
ent lithologies. 

A second color composite for the test site is presented in 
Figure 67. The same three components and color assignments as the last figure 
were used, except that the uncomplemented visible image was used for the green 
component (spectrally dark areas are dark, bright areas have a high green 
content). The accompanying interpretation map (Figure 68) is very similar to 
Figure 67, with a few exceptions. Differences are now discernible in the 
Death Valley salt pan - some areas are green, others are cyan. The inter- 
mountain valleys and playas also display differences not previously discern- 
ible. These differences are attributable to changes in material (salt vs. 
clay vs. sand) and to moisture content. 

C. Plsgah 

1. Geologic Description 

The Plsgah test site is located in the Mojave Desert region of 
southeastern California (Figure 48). Topography is variable: flat alluvial 

fans and dry lake beds surround rugged mountains with relief in excess of 
700 m. Vegetation Is either sparse or absent owing to the minimal annual 
prec ipltation. 

Rocks exposed In the area are igneous intrusives and extru- 
sives, and sediments derived from these sources. The intrusive rocks include 
quartz monzonlte, blotlte quartz monzonite, felsite, and mafic and felsic 
dikes. Extrusive rocks Include basalt, andesite, dacite porphyry, and pyro- 
clastic assemblages of variable lithology. Sedimentary rocks are dominated by 
fanglomerates and fan gravels derived from adjacent mountains and deposltrd as 



Figurt' (i7 . IK;MM color composite of Death Valley. Day IR, night IR, aiul 


visible Images displayed as hlvie, red, and green, respective- 
ly. Day and night IR Images were complemented. 
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DEATH VALLEY, CALIFORNIA: HCMM 



INTERPRETATION OF HCMM COLOR COMPOSITE: BLUE - DAY IR 
GREEN - DAY VISIBLE (4), RED - NIGHT IR (-) 


Figure 68. Interpretation map of Figure 67, 
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alluvial fans. Alkaline clay occurring In the dry lake beds and windblown 
sand deposits form surflclal sediments of non-volcanlc origin. 

Structures In the area are dominated by northwest trending 
right-lateral strike-slip faults with some normal displacement, and minor 
normal faults of various trends (Geologic Map, Figure 69). 

2. HCMM Images 

The HCMM day IR Image Is shown as Figure 70. The cold areas 
(dark) correspond to Lavlc Lake playa, and outcrops of pyroclastlcs, andesite, 
some basalts, daclte, blotlte granite, and granite. The area of aa at the 
Plsgah flow Is colder than the pahoehoe, though the flow Is not separable from 
the Surrounding alluvial fans. The warmest areas are those corresponding to 
older fan deposits (north central part of the Image). 

The night IR Image is shown as Figure 71. The warmest areas 
generally outline bedrock outcrops; alluvial areas are colder; Lavlc Lake 
playa Is quite cold. Very ; r 'le lithologic discrimination Is possible based 
on this Image. 

The visible Image (Figure 72) allows separation of several 
units. The darkest areas are basalt outcrops, playas are the lightest, meilum 
grays correspond to the rest of the Intrusive rocks. Again, very little 
discrimination Is possible, though there appears to be more information 
present here than on the night IR Image. 

In order to examine the Interrelationship of the three HCMM 
data types, several color conqsosite Images were produced by using various 
combinations of data types as components of a color additive triplet. One 
such Image is shown In Figure 73 (left). This Image Is a combination of the 
day IR, night IR, and visible, displayed as blue, red, and green respectively. 
All three conqjonents were complemented, so areas that are cold or have low 
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I'lfture 70. HCMM day IK Imago, of PLsgah aroa. Imago size is 
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HCMM nlHht IK Image of Pisgah area. Image size i.s 35 x 40 km. 
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Flmire 72. HCMM visible image of Pisgah area. Image size is 35 x 40 km. 
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alhudq appear aiaf strongly colored In the Image. For comparison, a landsat 
scene of the same area was computer processed; o£ several enhancements 
produced, the best, a principal components transformed Image (Blodget et al., 
197F), Is shown In Figure 73 (right). Interpretation maps of each are 
presented in Figure 74. 

Hie HCMM Image presents a wide variety of different colors for 
bedrock areas, whereas most of these areas are similarly colored on the 
landsat image. Many of the rocks a\ ?. quite flat spectrally in the landsat 
wavelength region, and so cannot be separated. These same rocks, however, 
appear to have different thermal properties, and so can be separated using 
HCMM data. 

On the HCMM Image, basalts are green. The aa flows at Pisgah 
and west of Pisgah are yellow, and easily separable. They are colder at night 
than the pahoehoe flows due to their lower thermal inertia. Granitic rocks 
are generally orange; the blotite quartz monzonites In the south are somewhat 
redder. The daclte west of Pisgah is the same orange color; composltionally 
this rock is similar to the granites. The sand dunes in the northwest come;;’ 
are red (cold at night, warm in the day, high albedo) as is La vie lake playa. 
The other extrusive rocks are orange-yellow in color. The varied colors of 
the alluvial areas reflect the variation in their source rock composition. 

On the Landsat Image the bedrock areas are easily distinguish- 
able from the alluvial areas. However, very few separations can be made 
between the various rock types. The major improvement over HCMM data Is the 
higher spatial resolution which permits mapping of smaller features. Very 
little other information can be interpreted from the Landsat Image. 
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3. Alrcrairt Images 

The greatly Increased spatial resolution provided by aircraft 
scanners is shovm in Figure 75, a thermal inertia image of the Plsgah flow and 
vicinity. An interpretation map produced from this image (Figure 76) and the 
description of the symbols (Table 20) are also presented. 

TABLE 20 

DESCRIPTION OF SYMBOLS FOR PlSGAH AIRCRAFT INTERPRETATION 


Unit 

Geologic 
Map Symbol 

Material 

1 

Qa 1 , Qco 

Alluvium 

2 

Qp 

Playa 

3 

Qb, Tb 

Basalt 

3a 


"Aa" 

4 

Ti 

Rhyolite intrusive 

5 

Gr 

Granite 

6 

Bm 

Blotite quartz monzonlte 

7 

Mvp 

Pyroclastics 

8 

Mva , Mvb 

Andesite, basalt 


Some thermal inertia differences associated with different rock 
types are depicted by varying gray tones in the image. The area of aa flows 
(#3a) is clearly separable from the pahoehoe flows (//3). The low thermal 
inertia area of lavic Lake playa In the south central part of the image Ip 
also readily apparent. A major characteristic of the terrain, discernible due 
to the higher resolution, is geomorphology and drainage texture. With the 
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Figure 75. Thermal inertia image of Pisgah from aircraft data. Imag..* size 
IS 18 X 20 km. 
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Figure 76, Interpretation map of Figure 7.,, 
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addition of this Information, several units can be separated. The older fans, 
for example, have a dissected appearance due to high drainage density. The 
rhyolite Intrusive (M) appears distinctive due to linear features (dikes) 
cutting It. Plow features are visible on the Plsgah flows. Alluvial fans 
appear fairly homogeneous with parallel to braided drainage patterns. This 
type of Information, vital to photointerpreters. Is not available from HCMM 
data. 

The north-northeast trending linear feature In the southwest 
part of the Image (indicated by arrows) is interpreted to be a fault, and does 
not appear on the detailed geologic map of the area. It is easiest to distin- 
guish this feature on the thermal Inertia Image. It seems to be expressed as 
a narrow valley or series of low areas covered with alluvium, which have a 
lower thermal inertia than the surrounding bedrock. Future field work will be 
necessary to determine the actual cause of this feature. 

D. San Rafael 

1. Geologic Description 

The San Rafael Swell is a large asymmetric anticline located In 
Emery County, Utah. (Figure 77;. Rocks ranging in age from Pennsylvanian to 
Cretaceous are exposed. The Swell is a breached structural feature: the inner 
part is underlain mainly by Permian and Tr lassie age rocks; the outer part is 
formed mainly of Tr lassie and Jurassic age rocks; younger strata of Jurassic 
and Cretaceous age ring the structure (Hawley, et al. , 1968). 

The oldest rocks, of Pennsylvanian and Permian age, are exposed 
in deep canyons and consist of the Elephant Canyon formation and the Hermosa 
formation. Younger strata in the Inner part are more widely exposed and 
consist of the Coconino Sandstone and Kalbab Limestone of Permian age, and the 
Moenkopl and Chlnle formations of Triassic age. In the outer part of the 
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Figure 77. Location map of San Rafael Swell, Utah, 
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Swell are found exposures of the Wingate Sandstone, Kayenta formation, and 
Navajo formation of Triassic and Jurassic age. Ringing the Swell are ex- 
posures of the Carmel formation, the Entrada Sandstone, the Curtis formation, 
the Summerville formation, and the Morrison formation of Jurassic age, the 
Cedar Mountain formation, the Dakota Sandstone, and the Mancos Shale of 
Cretaceous age (Figure 78). 

Vegetation throughout the San Rafael Swell Is generally sparse, 
ranging from nearly 0% to 20% atop limestone and sandstone benches. At higher 
elevations. In mountains surrounding the Swell, substantial plnon- juniper 
forests are found with 80-100% vegetative cover. 

2. HCMM Image 

The HCMM day IR Image for the San Rafael Swell Is shown In 
Figure 79. An Interpretive map derived from this Image, and explanation are 
presented In Figures 80 and 6a. 

Overall, this Image allows separation of only a few geologic 
units. The most consistently distinct unit is the Entrada Sandstone, a red, 
fine grained sandstone which weathers to produce sandy/sllty flat-lying areas. 

The Nava jo/Kayenta/Wingate Sandstones are usually separable as 
a group. They are. In general, cllf f-forralng, massive, light-colored sand- 
stones, distinct from the overlying Carmel Formation and underlying Ghlnle 
Formation. 

Separation of most of the other units was difficult or impos- 
sible due to the limited spatial resolution of the HCMM data and the relative- 
ly narrow outcrop width of many of the geologic units. 

Tlie night IR HCMM Image (Figure 82) provided some improvement 
for separating various units (Figure 83). The massive light sandstones appear 
warm In the Image; the Entrada appears cold. Contrast between other, adjacent 
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Figure 78. 


Stratigraphic column of rocks exposed in San Rafael Swell. 
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i-'igure 79. IICMM day IH imago of San Rafaol Swoll. Imago size is 140 x 14 < .• n. 
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Figure 80 . Interpretation map for Figure 79 
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Ktratlgraphic units allows tham to be separated. In general, warmer areas are 
those underlain by massive, light sandstones (high thermal inertia), and dark, 
cold areas are underlain by silty material. The effect of elevation on 
temperature is indicated by the cold ridges in the northwest corner of the 
image. 

The thermal inertia image (Figure 84) generally provided the 
same amount of information (Figure 85) as the day or night IR images. Bright 
areas (high thermal inertia) are underlain by the massive sandstones. In the 
northwest and northeast corners of the image, the high thermal inertia values 
are due mainly to high elevations (no corrections for elevation-dependent 
temperature variations were made) and cold nighttime temperatures. Low 
thermal inertia areas are underlain by unconsolidated material such as the 
sand dunes in the southeast part of the image, and the Entrada and Carmel 
Formation benches. 

Comparison of the day and night IR data indicated that there 
was information on each not shown on the other. In order to combine the data, 
simple color additive composites were produced. The best one is s^hown in 
Figure 86. This composite has as its components day IR, night IR, and visible 
displayed in red, blue, and green respectively. All three components were 
complemented: areas that are either spectrally dark, or cold are displayed as 
brightly colored on the image. 

Generally, more units could be separated using the color 
composite than any of the single images previously described (Figure 87). 

Along the west side of the Swell, the Tununk shale, Dakota Sandstone, and 
Cedar Mountain, Brushy Basin, Entrada, Carmel, Navajo and Chinle formations 
were separable. These units have appropriately large outcrop patterns, vary- 
ing albedos, and varying thermal properties. The Entrada Sandstone, for 
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ligure 84. HCMM thermal inertia image of San Rafael Swell. 
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Figure 85. Interpretation map for Figure 84. 
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f'igure 86. HCMM color composite of San Rafael Swell, Day IR, night IR, and 

visible images displayed as red, blue, and green, respectively. All 
three components were complemented. 
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example, ts dark blue, which indicates It has low albedo, is warm in the day, 
and relatively cold at night implying a low thermal inertia. The Navajo sand- 
stone is orange on the image, which indicates it has high albedo, is cold in 
the day, and warm at night. This indicates a high thermal inertia. Similar 
conclusions can be made about the other units separable on the image. The 
appearance on the image of the units is consistent with what is known about 
their characteristics. 

The color compositing technique presented here appears to be a 
simple but effective method of displaying the various HCMM data types acquired 
by the satellite. Tlie use of color Improves the analyst's ability to recog- 
nlase differences between areas, and allows the simultaneous display of multi- 
variate data. 

For comparison purposes, a landsat-Seasat composite scene 
covering the same area was obtained and interpreted (Figure 88). This image 
is a combination of coregistered Landsat and Seasat data. Processing consist- 
ed of performing a principal components transformation of the Landsat bands , 
then using the Seasat radar data to modulate the intensity of each of the 
three components for i:he color composite. In this manner the higher resolu- 
. ion of the Seasat data is merged with the Landsat spectral data, with some 
contribution from Seasat backscatter response folded into the color picture as 
intensity (brightness). This product was found to be the most effective 
combination of several examined. 

The interpretation map (Figure 89) shov;s that significantly 
more lithologic separations can be made using Landsat/Seasat data compared to 
HCMM data. Tills is due both to the improved spatial resolution, and to 
spectral bands of Landsat and the Seasat backscatter response. Contrary to 
the results from some of the other test sites, no additional information was 
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Figure 88. Landsat/Seasat composite of San Rafael Swell, Utah. 
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Figure 89. Interpretation map for Figure 88. 
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obtained from the HCMM data that was not already available In the Landsat/ 
Seasat data. 

The spectral characteristics of the sedimentary .ocks exposed 
In the San Rfifael Swell are sufficiently distinctive to allow their separation 
In the Landsat wavelength region. In addition, textural Information Inter- 
pretable from the Im/Ages, alec contributes to the separability of the units. 
However, this critical characteristic cannot be extracted or Inferred from the 
HCMM data due to the poorer spatial resolution. 
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VXI. SUMMARY AND CONCUUSCONS 


Our attum|>ts to determine quantitatively accurate thermal inertia values 
from HCNM, digital data have met with only partial success owing to a variety 
of factors. Calculations of the atmospheric radiative transfer effects showed 
that the radiance measured at the satellite was equivalent to a surface tem- 
perature of from f3 to -ll^C different from the actual ground temperature. 
These effects were strongly dependent on water vapor content of the atmos- 
phere, surface elevation and surface temperature. These effects could easily 
produce an error in AT of 6*^0, and perturbation studies with our thermal model 
showed that even a error In Al' could cause an error In thermal inertia of 
0.01 cal . Because of the delay In receipt of HGMM data and the 

continued confusion over the calibration of the HCMM sensor, we were forced to 
leave Incorporation of these importa, t atmospheric effects to a future study. 
The effects of atmospheric transmission and radiance, instrument calibration, 
and varying elevations are tractable in future missions and some success 
should still be attainable upon ' rther work with HCMM data. However, a 
solution to the problem of varying meteorology across a scene (as it affects 
the surface heat balance, not as required for atmospheric radiance correc- 
tions) may prove to be prohibitively costly and time consuming. 

As an adjunct to the satellite and aircraft studies, we built and tested 
a thermal inertia meter (TIM) In order to obtain in situ thermal inertia 
values of rock and soil types. The results of laboratory tests indicated that 
precise thermal Inertia values were not possible to achieve principally 
because of the difficulty in providing a constant heat flux it the surface of 
the standard and target materials. Nonetheless, field tests in Nevada and 
California desert regions showed that the TIM was useful in a relative sense. 
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As lortg as wind speeds were low, to minimize convective heat loss, TIM mea- 
surements were accurate enough to be used to discriminate many dllferent types 
of geologic materials from one another, e»g< chert and barite* 

Spatial and temporal soil moisture effects are evident In the therroal 
Inertia images constructed from both the aircraft and HCMM data. Parts of 
Death Valley's floor are subject to seasonal flooding and the lower half of 
the Valley, around Badwater Basin, Is nearly always saturated* On the thermal 
Inertia image, these moist areas are bright corresponding to regions of rela- 
tively high thermal Inertia values* Similar patterns are evident on the dry 
lakes In the Walker l^aine area of Nevada, and some evidence of seasonal change 
can be seen when comparing Image.^ at different times of the year. The details 
of seasonal moisture changes are important for future studies but Vtave not 
been treated here in any depth. This requires additional field work and the 
uoe of the soli moisture version of the JPL thermal model. This model has not 
been used to analyze the HCMM data for this study, but It Is available now for 
a future effort* 

An Important positive result that has come out of this study Is the 
recognition that In most Instances apparent thermal Inertia Is an excellent 
qualitative representation of true thermal Inertia* For many applications 
this will considerably reduce data processing requirements, with only minimal 
loss of Information. Apparent thermal Inertia should prove to be completely 
adequate for those geological applications where the desire Is to use the data 
as part of a photo-geologic Interpretation problem such as to discriminate 
between units or to examine structural or textural relationships* Little loss 
of Information should be noticeable when using apparent thermal Inertia in 
combination with other data sets, such as Landsat data. In multivariate 



analysis alf?ottthms , as long «s the analysis Is confined to small, local 
studies* 

The apatlal resolution of HCHM Is considerably less than desirable for 
tK‘ type of geologic Investigations which we hove undertaken. Nevertheless, 
considerable success has been achieved in demonstrating the utility of day- 
night thermal data for geologic applications. For the four sites examined, 
Only limited geological information was Interpretable directly from Images of 
individual HGMM data types; l.e., day visible, day IR, night IR, or thermal 
inertia. Various display enhancetnent techniques were applied to these data to 
try to Increase the information displayed, with little success. A method of 
combining HCMM data types to create color composites was developed which: a) 

was computationally simple; b) allowed simultaneous display of several vari- 
ables; c) effectively displayed the Interrelationship of the HCMM variables; 
d) provided a useful, interpretable product displaying significantly more 
lichologlc information. 

At the Walker Lane site, the color composite of HCMM data provided some 
lithologic separations. Playas and alluvial areos were generally distinguish- 
able; the alteration zones at Goldfield and Cuprite had similar appearances to 
the playas; basalts and other volcanic rocks were mappable, though there was a 
significant proportion of the scene where these rocks were not separable; 
large areas were found where no useful lithologic distinctions could be made. 
The high resolution aircraft data allowed identification of the slllclfled 
altered rocks at both Goldfield and Cuprite due to th^ir high thermal inertia. 
A combination of the thermal inertia data with visible and near Infrared 
multispectral data was produced for Cuprite. The thermal Inertia data 
Improved the separation of sub-units of unaltered, spectrally dark and flat 
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volcanic rocka* No other addlclonol InforiMtion waa dlacernible due to the 
addition of thermal data. 

At Death Valley > HCMM thenwil Inertia imagea allowed delineation of oreaa 
underlain by dolomite, llmeatone, quartzite, and granite (all high thermal 
inertia materlnla). The boundary between bedrock and alluvlvim waa clearly 
dlac'rnlble and eaaler to recognize than from Landaat Images. Examination of 
coll r compoBltea of HGMM data allowed separation of various lithologies and 
alluvial units. Playas, saltpans, sand, sediments, and volcanic rocks were. 

In places, separable. 

At the Plsgah site, the Individual HCMM data types provided limited 
separation of lithologic units. The color composite produced using theoe data 
types provided a significant amount of information. The basalt flows of 
Plsgah could be separated Into aa and pahoehoe units; other basalts were 
distinctive in color; granitic rocks and ml nera logically similar volcanic 
rocks were separable, and had tilmllar appearance on the Images; boundaries of 
various types of alluvium were distinct. A comparison with a computer pro- 
cessed Landsat scene Indicated that despite the higher spatial resolution of 
Landaat data, HCMl data provided more lithological Information. This Is 
thought to be due to the featureless spectral properties of the materials in 
the landsat wavelength region versus the variable thermal properties In the 
HCMM wavelength region. Examination of high resolution aircraft thermal data 
revealed the presence of a previously unrecognized linear feature, interpreted 
to be a fault. It was most easily discernible on the thermal inertia image. 

Finally, at San Rafael Swell a comparison of HCMM data (both Individual 
bands and color composites) with a Landsat/Seasat composite Indicated that no 
additional lithologic Information was provided by the HCMM data. The spectral 
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characteristics of the units and the Improved spatial resolution of Landsat 
and Seasat allowed discrimination of almost all of the mapped geologic units. 

First attempts at digitally combining HCMM data with Landsat multi- 
spectral data at the various sites show promise for Improved rock discrimin- 
ation. In a future study we Intend to continue to evaluate the Improvement In 
geologic remote sensing made possib'ie by this extended multivariate data set. 
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